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ABSTRACT 

We present a systematic study of the shape of the dust attenuation curve in star-forming 
galaxies from the far-ultraviolet to the near infrared (~ 0.15 — 2/xm), as a function of specific 
star formation rate (tps) and axis ratio (b/a), for galaxies with and without a significant bulge. 
Our sample comprises 23,000 (15,000) galaxies with a median redshift of 0.07, with photo- 
metric entries in the SDSS, UKIDSS-LAS (and GALEX-AIS) survey catalogues and emis- 
sion line measurements from the SDSS spectroscopic survey. We develop a new pair-matching 
technique to isolate the dust attenuation curves from the stellar continuum emission. The main 
results are: (i) the slope of the attenuation curve in the optical varies weakly with ips> strongly 
with b/a, and is significantly steeper than the Milky Way extinction law in bulge-dominated 
galaxies; (ii) the NIR slope is constant, and matches the slope of the Milky Way extinction 
law; (iii) the UV has a slope change consistent with a dust bump at 2175A which is evident in 
all samples and varies strongly in strength with b/a in the bulge-dominated sample; (iv) there 
is a strong increase in emission line-to-continuum dust attenuation (Tvjine/Tv.cont) with both 
decreasing ips an d increasing b/a; (v) radial gradients in dust attenuation increase strongly 
with increasing ips> an d the presence of a bulge does not alter the strength of the gradients. 
These results are consistent with the picture in which young stars are surrounded by dense 
'birth clouds' with low covering factor which disperse on timescales of ~ 10 7 years and the 
diffuse interstellar dust is distributed in a centrally concentrated disk with a smaller scale- 
height than the older stars that contribute the majority of the red and NIR light. Within this 
model, the path length of diffuse dust, but not of birth-cloud dust, increases with increasing 
inclination and the apparent optical attenuation curve is steepened by the differential effect of 
larger dust opacity towards younger stars than towards older stars. Additionally, our findings 
suggest that: (i) galaxies with higher star formation rates per unit stellar mass have a higher 
fraction of diffuse dust, which is more centrally concentrated; (ii) the observed strength of 
the 2 175 A dust feature is affected predominantly by global geometry; (iii) only highly in- 
clined disks are optically thick. We provide new empirically-derived attenuation curves for 
correcting the light from star-forming galaxies for dust attenuation. 
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1 INTRODUCTION 

Dust is an ubiquitous component of the baryonic Universe, po- 
tentially condensing out of the interstellar medium as soon as 
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the first ge neration of stars has expired and provided the requi- 
site metals (Meurer et al.l ll997l: iHughes et alii 19981 : iBertoldi etafl 
120031 : iDraind 12009). The presence of dust in galaxies hampers 
our measurements of intrinsic galaxy properties, including stellar 
masses and star formation rate (SFR), when observations are lim- 
ited in wavelength coverage or availability of spectroscopic data 
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dCalzetti et al.l2007h . It also impacts on our measurements of global 
population statistics, such as luminosity functions and colour- 
magnitude diagrams. The strong systematic variation of dust con- 
tent with stellar mass, SFR and morphology can lead to significant 
biases in results if not correctly accounted fo r ( Brinchmann et al 



20041; IDriver et al.l2007l : lda Cunha et al.l2010l : lGilbank et al.l201 



Mailer et al. 2009). Dust expelled from galaxies and re siding in 
galaxy halos dWild & Heweta20o3 : lHewett & WildteOOTl) can even 
have implications f or the estimation of cosmological parameters 
dMenard et ai1l20ld) . 

For several decades research into dust in galaxies focused 
around the question of whether spiral disks are optically thick or 
thin. By studying large samples of galaxies with varying incli- 
nations to the line of sight, in combination with models for the 
dust distribution in galaxies, NIR photometry, or rotation curves, 
the amount of attenuation suffered b y starlight before it reaches 
the observer can be constrained (e.g . [Giovanelli & Havnesll2002l ; 
iDriver et al.l2007l ; iMaller et alj2009tlYip et alfeOlOh . The consen- 
sus is that the stellar continua in galaxy disks suffer 0-1 magnitude 
of attenuation in B (or SDSS (/-band), dependent on inclination. 
Galaxy bulges may well suffer greater a ttenuation due to radial gra- 
dients in dust disks. [Driver e t al. (2007) find inclination dependent 
attenuations in bulges of 0.8-2.6 magnitudes in the B-band. 

While the extinctiorQ of starlight by dust can be measured 
with relative ease along lines of sight to stars in the Milky Way 
and closest galaxies (SMC and LMC), the integration of light over 
the entire stellar population of a distant galaxy introduces con- 
siderable complications. Firstly, the dust is not distributed uni- 
formly throughout the galaxy: there are dense birthclouds around 
the youngest stars, th ere is a diffuse component exhibiting strong 
radial gradients (e.g . [Persson & Helou 1987; Pe letier et al.ll 199a ; 
iBoissier eta l. 2004; Munoz-Mateos et al. 2009), and there is some 
evidence that diffuse interstellar dust is distributed t hroughout the 
disk w ith a smaller scaleheight than that of the stars dxilouris et al.l 
1999). Secondly, the dust does not approximate a uniform screen 
in front of the stars: dust and stars are mixed, potentially causing 
an additional scattering component which can affect the shape of 
the atte nuation curve r elative to the underlying extinction curve. 
Indeed, iBoauien et all (2009) suggest that variation in extinction 
laws and dust-star geometry may cause the scatter in the dust atten- 
uation vs. re ddening (e.g. IRX -/3) relation observed in star-forming 
galaxies, and Bu at et alTd201ll) find tentative evidence for a steeper 
dust attenuation curve i n ordinary star-formi ng galaxies than found 
in starburst galaxies bv lCalzetti etal] d2000h . 

On the other hand, the sensitivity of the shape of the dust 
attenuation curve to dust grain properties and dust-star geome- 
try allows us to extract information on the properties and distri- 
bution of dust in galaxies. For example, a flatter blue/UV atten- 
uation curve is generally expected in the case of spatially mixed 
dust and stars, as the least extincted stars contribute a greater frac- 
tion of total emergent blue light than red light (e.g. iGordon et al.l 
1 19971 ; ICharlot & Fall 2000). Dust features, the strongest one be- 
ing at 2 175 A, can provide further constraints. Models suggest that 
the 2175A feature should vary in strength in galaxy spectral en- 
ergy distributions (SE Ds) depending upon the properties and ge - 
ometry of the dust dGranato et all 12000); Iwitt & Gordonl |2000|) . 



Observationally, there is some evidence for a 2175A abso rption 
feature in ordinary sta r-forming galaxie s dConrov et al.ll20 100 and 
high-redshift galaxies (Noll et al. 2009), but clear absence in local 
starbur st galaxie s (Calzetti et al. 1994) and the Small Magellanic 
Cloud dPeill992l ; lGordon et al.l2003l) . 

In this paper we measure empirically the optical-NIR and UV- 
to-NIR attenuation curves, using a sample of ~23,000 and ~15,000 
star-forming galaxies respectively, as a function of radius, specific 
star formation rate (ips= SFR/M*), inclination angle (b/a), for 
galaxies both with and without massive bulges. We achieve this by 
carefully pair-matching galaxies to remove unwanted stellar popu- 
lation signal. Spectroscopic information provides a measure of rel- 
ative dust content via the Balmer emission lines, allowing us to 
combine galaxy pairs in such a way as to reveal an average attenu- 
ation curve with very high SNR. 

The development of this new, robust, statistical method to 
measure dust attenuation curves in galaxies with a range of dif- 
ferent underlying stellar populations, allows us for the first time to 
observe trends in the shape of the attenuation curve as a function 
of galaxy properties. Deriving attenuation curves using different 
sized apertures, enables us to measure the strong radial gradients in 
dust opacity with tps- We compare dust opacity suffered by emis- 
sion lines to that suffered by the optical continuum, finding strong 
trends with tps and axis ratio. Finally, we study the strength of the 
2 175 A dust bump which we find to depend on inclination and the 
presence of a bulge. All these empirical results relate to the spatial 
distribution of dust in galaxies, relative to stars of different ages. 

In Section [2] we present our methodology, along with a brief 
summary of the formalism used to describe dust attenuation in 
galaxies. In Section [3] we describe the datasets used and provide 
details of the photometric information extracted from each one. In 
Section [4] we present the attenuation curves, the shape of which 
we characterise by measuring the slope in small wavelength ranges 
in Section [5] In this section we also present the radial gradients, 
a comparison of attenuation in the emission lines and stellar con- 
tinua, and the strength of the 2175A dust bump. In Section[6]we 
provide an empirically derived fitting formula to correct the spec- 
tral energy distributions of a wide range of star-forming galaxies 
for dust attenuation. Finally in Section [7] we present a qualitative 
model for the local and global geometry of dust in galaxies, which 
is consistent with the majority of our results. 

Throughout the paper we compare our results to the Milky 
Way (MW), large and small Magellanic cloud extinction curves 
(LMC and SMC), and the local starburst attenuation curve of 



Calzetti et alj <200CL COO). We adopt the MW extinction curve of 
Cardelli et al.l d 19891) . which is defined over the wavel ength range 
0.125 < \/fim< 3.5, with the correction given by lO'Donnelll 
dl994l) in the optical region 0.3 < X/ um< 0.9. We adopt the tab- 
ulated LMC and SMC-bar curves of bordon et alj d2003l) . spline 
interpolating onto our wavelength points. 

Absolute quantities such as stellar mass and star formation rate 
are required to match pairs of galaxies, and for these we assume the 
standard ACDM cosmology with Qyi = 0.3, J1a = 0.7 and h = 
0.7. The attenuation curves are based on flux ratios, and therefore 
independent of cosmology, magnitude types and filter zero-points. 



We follow the standard terminology (e.g. ICalzettil llooih : extinction = 
absorption + scattering out of the line-of-sight; attenuation = absorption 
+ scattering in to and out of the line-of-sight caused by local and global 
geometric effects. 



2 METHODOLOGY 

The basic principle of our method is the same as that used to 
measure the dust extinction along lines-of-sight through the Milky 
Way and Magellanic Clouds, where the spectral energy distribu- 
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tions (SED) of pairs of stars of similar spectral type are compared. 
Here, we calculate the ratio of pairs of galaxy SEDs, where one 
galaxy is more dusty than the other, combining many pairs to ob- 
tain a high SNR measurement of the attenuation curve. Just as with 
sight-lines to stars in the Milky Way or local galaxies, the method 
is only valid where no significant bias exists between the spectral 
types of the more dusty and less dusty galaxies. To ensure this, we 
match galaxy pairs in gas-phase metallicity, specific star formation 
rate (ips, corrected for dust attenuation using the ratio of Ha to H/9 
luminosity), axis ratio and redshift. 

To ensure that we always divide the SED of the dustier mem- 
ber of the pair with that of the less dusty member, we use the ratio 
of Ha to H/3 luminosity as a measure of relative dust content. Our 
method is independent of, and insensitive to, the absolute dust con- 
tent of the galaxies. In the first part of the paper we calculate the 
average attenuation curve for samples split by specific star forma- 
tion rate (ips), axis ratio (b/a) and stellar surface mass density (n*), 
in circular apertures of increasing physical size. We use the result- 
ing attenuation curves to study the changing shape and amplitude 
of the attenuation curve with each of these properties. In the second 
part of the paper we provide an empirically derived prescription for 
the attenuation of starlight by dust in external galaxies as a function 
of their ips ,b/a and yt* . 

2.1 Background to dust attenuation laws 

The many different formalisms used to describe dust curves neces- 
sitate a brief overview of the subject to allow orientation of our 
work relative to other papers in the literature. 

Light emitted by stars in external galaxies is affected by dust 
in three different ways: (i) absorption; (ii) scattering out of the line- 
of-sight to the observer; (iii) scattering into the line-of-sight to the 
observer. The first two combined are usually termed "extinction". 
In this paper we use the term "attenuation" to mean the average 
loss of light after integration over all lines-of-sight to light sources 
throughout the entire galaxy, encompassing all effects arising from 
the local and global geometric configuration of the dust and starfl 
We use the term "effective optical depth" to describe the final ob- 
served optical depth, which is related to the attenuation through 
A\ — 1.086ta. Because blue light is scattered and absorbed more 
than red light by dust grains, an overall "reddening" of the light is 
observed. In this section, we introduce the basic formalism used to 
describe the ef fects of dust on the integrated light from galaxies. 
ICalzettil ( feOOll) gives an extensive review of the literature up until 
a decade ago, although note that the terminology differs in places 
from that adopted here. 

The attenuation A\, in magnitudes, at a given wavelength A is 
given by: 

A x = -2.51og 10 (^ (1) 

= -2.51og 10 [exp(-TA)] (2) 
= 1.086ta (3) 

where is the observed luminosity, I\ is the intrinsic (emitted) 
luminosity of the source, and Ta is the effective optical depth of the 
dust. Ultimately it is this quantity that we wish to know in order to 
correct the SEDs of galaxies for dust attenuation. 

2 Some other papers refer to "single-line-of-sight" attenuation, which in- 
cludes effects (i)-(iii) but excludes global geometric effects. 



Observationally, the measurement of the overall amplitude of 
extinction and reddening of the light are usually separate tasks. It is 
therefore practical to separate the attenuation into two components, 
the amplitude (usually specified in the V-band, ty) and the shape 
of the attenuation curve (<2a0 

Ta = r v Qx (4) 

where Qv = 1. In the case of extinction, Q\ is a function of 
the distribution of grain composition, size and shape. When study- 
ing the integrated light from external galaxies Q\ is additionally a 
function of the local geometry of the dust and stars, and the global 
geometry of the galaxy (bulge and disk, and inclination). 

Extinction curves have been measured along sight lines to- 
ward stars in the MW, LMC and SMC. While the average curves 
of the three galaxies are apparently entirely distinct, the individ- 
ual sight-lines in all three galaxies show a continuum of properties 
from steep slopes with weak 2 175 A dust bumps like the average 
SMC- bar curve, to shallow slopes with strong bumps like the MW 
curve l lGordon et al.ll2003h . Unfortunately, current data is not good 
enough to determine robustly the "hidden parameter(s)" responsi- 
ble for this variation, although grain destruction by local star for- 
mation is often assumed. 

The attenuation suffered by light from stars in nearby star- 
burst galaxies, where the unattenuated stellar continuum varies rel- 
atively little between galaxies, has been measured empirically from 
the ultra-violet (UV) to near infra-red (NIR) i n a series of papers 
by D . Calzetti and co l laborators (Kinney et alJI 19 94; Calzet tiet al.l 
1 1994 ICalzettill 19973 : iBoker et al.ll 19991: ICalzetti et alj|200d) . The 
attenuation curve was found to be "greyer" (i.e. flatter) in the opti- 
cal than the Milky Way extinction curve, qualitatively as expected 
from the geometrical effects of mixed dust and stars. Additionally, 
no evidence was found for a strong "2175A" dust feature which is 
observed in the majority of ex tinction curves measured along lines- 
of-sight in the MW (but see iNoll et alj|2009l : IConrov et al j|201Ct 
and Section |53t . 

2.2 Measuring dust attenuation curves from pair-matched 
galaxy samples 

An observed galaxy SED can be described as an intrinsic SED of 
unit stellar mass (1%), multiplied by factors to account for the to- 
tal stellar mass, cosmological distance of the galaxy and internal 
dust attenuation. Ideally the two factors of normalisation should be 
estimated from a region of the SED that is not attenuated by dust. 
In this paper, the longest wavelength point available to us is the 
K-band at 2.2/im. We therefore normalise each galaxy by its K- 
band flux and retain a small correction factor throughout this work 
which allows for the small amount of dust attenuation in the K-band 
(Atk) during the fitting procedures. The normalised SED (Fx) can 
be written as: 

F A = ^exp(-T A ) (5) 

JK 

where A accounts for stellar mass and cosmological distance, 
exp(— Ta) for internal dust attenuation, and /« is the observed K- 
band flux. 

Directly from Eqn. [5] the ratio of a pair of normalised galaxy 
SEDs (labelled 1 and 2), where both members of the pair have the 

3 Conversions between different parameterisations of attenuation curve 
shapes are : E{ ^_ v) E{B A ~ V) =^ = ^ = Qx 
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same intrinsic SED but different dust contents and absolute normal- 
isations, can be written after some algebra as: 

= exp(-AT A + At k ) (6) 

where Ata = Ta,i — Ta,2 and Atk accounts for the small, un- 
known extinction in the K-band. The geometric mearQ of a sam- 
ple of flux ratios, where in each pair galaxy F 1 is more dusty than 
galaxy F2 , results in a dust attenuation curve normalised to zero at 
K which we denote 7a to distinguish it from the true attenuation: 

= - (At v > Qa + (At k ) (8) 

where n is the number of pairs of galaxies in the sum, (Ata) is 
the arithmetic mean of the difference in dust attenuation between 
galaxies 1 and 2, and we have expanded Ata into the amplitude of 
attenuation at 5500A (Atv) and the shape of the attenuation curve 
(Qa), as described above. From now on we will drop the brackets 
and write simply At as the effective optical depth averaged over 
multiple galaxy pairs. 

Some points to note about the method are: 

• Because we are working with broadband photometry, strictly 
Eqn. [5] should include an integration over the filter response func- 
tion (R\). In the case that Qa varies with A over the wavelength 
range of the filter, it can be shown that we measure J Q\R\d\ 
rather than Qa, even in the case that IJ varies with A. Due to the 
generally weak variation of attenuation curves with A, this distinc- 
tion is only important in the UV (Section |53) and we have therefore 
omitted the integrand for clarity. 

• Eqn.[8]holds only in the case where Qa is the same for both 
members of the pair i.e. when there is no trend between Qa and 
dust content. See Appendix lAl and Figure lATl for the result of tests 
showing that this is the case in our samples. 

• Qa may be different for different pairs in the sum, in this case 
the measured Qa will be the Atv weighted mean Qa of the sample. 

• The method does not require that the intrinsic stellar popula- 
tions are exactly the same in each galaxy pair, any residual stellar 
population remaining after division of the two SEDs will only add 
noise to the final measurement. However, it does require that there 
is no systematic correlation between the residual and the dust con- 
tent of the galaxy. See Appendix [A] and Figure lA2l for an example 
of how such a correlation can occur using photometry to estimate 
ipS rather than the Ha emission line. 

• The method is almost, but not quite, the same as that employed 
by Calzetti et al. for starburst galaxies: 

- They use a single "dust-free" composite spectrum with 
which to compare dustier composite spectra. Because we use a 
whole range of types of galaxies, with many different shapes of 
SEDs, we first define strict galaxy pairs, the SEDs of which we 
divide before combining galaxies into samples. This approach 
makes our method more robust, as each pair is treated individu- 
ally meaning outliers can be easily identified and removed. Ad- 
ditionally, errors can be estimated on the resulting attenuation 
curves. 



4 The geometric mean of flux ratios is preferred for dust curves , as this is 
equivalent to the arithmetic mean of the log of flux ratios, which is directly 
related to optical depth. 



- Each of our binned samples of galaxies contains 14- 
40 times more gala xies than the original starburst sample of 
ICalzettietai] ri994) which contained 39 galaxies. 

- We define the attenuation curve in terms of the shape of 
the curve (Qa), and the amplitude of attenuation in the stellar 
continuum at V (Ty). This differs in formalism to that of Calzetti 
et al. who directly incorporate the amplitude of dust attenuation 
in the emission lines into the attenuation curve for the stellar 
continuum. 

2.3 Pair-matching procedure 

In order to achieve the requirement for a close match in intrinsic 
stellar population, we identify pairs of galaxies that match in metal- 
licity, ips, b/a and redshift. The restrictions on metallicity and ips 
are to ensure that both galaxies have similar intrinsic SEDs. The re- 
striction on b/a also minimises the effect of seeing different stellar 
populations in inclined vs. face-on galaxies simply due to changes 
in transparency. 

In order to identify galaxies which are more- or less-dusty we 
use the ratio of Ha to H/3 luminosity measured from the SDSS 
spectra. At this stage, the sole purpose of using this line ratio as 
a measure of dust content is to ensure that, in as many pairs as 
possible, the denominator SED is more dusty than the numerator in 
Eqn. [7] Only in this case will we measure an overall attenuation. 

Because the Ha and H/3 emission lines are measured within a 
small fixed angular aperture (see Section |3~T) , and radial gradients 
in dust exist, both galaxies in the pair must have the same redshift 
to minimise the numbers of pairs in which the denominator SED is 
in fact less dusty than the numerator simply because the lines are 
measured from different regions of each member of the pair. The 
restriction on redshift also means that galaxy pairs do not need to 
be corrected for any relative shifting of band-passes with redshift, 
as the photometric bands probe the same rest-frame wavelength of 
the SED in each member of the pair. 

The procedure is as follows. The galaxies are ordered by in- 
creasing ratio of Ha to H/3 luminosity, and starting from the least 
dusty galaxy in the sample, the dustiest suitable partner is found us- 
ing the criteria mentioned above. We require that there is a signif- 
icant difference between the dust content of galaxy-pairs, as mea- 
sured by the difference in emission line optical depth. This max- 
imises the signal-to-noise of the final attenuation curves, by ex- 
cluding pairs which contribute only noise and no signal. 

To summarise the pair-matching criteria: 

• |AlogO/H| < 0.05 

• I A log tp s \ < 0.1 

• |A(6/a)| < 0.1 

• \Az\ < 0.01 

• AT V ,|ine> 0.2 

Details of how the derived quantities are measured are given in Sec- 
tion [3] If a suitable pair is not found, the initial galaxy is removed 
from the sample, and the algorithm moves onto the next galaxy. 
Once a pair has been identified, both members are removed from 
the sample, thus each galaxy can only contribute once to the atten- 
uation curve. 

Once all possible pairs have been identified within a particular 
sample, the galaxy SEDs are normalised by their K-band flux and 
the flux ratios calculated in each band for each pair. These flux ra- 
tios are combined using a geometric mean (Eqn. [7}. A preliminary 
stack is made, outlying pairs with flux ratios which are in the upper 
or lower 1% of the distribution in any band are removed from the 
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sample, and the attenuation curves recalculated using the cleaned 
sample. 

We verified that the precise details of the pair-matching cri- 
teria given above are unimportant to the final results by varying 
the widths of the bins and repeating our analysis. Tests showed 
that the best SNR was achieved when the difference between the 
dustiness of each pair of galaxies is maximised. We verified that 
selecting pairs to have approximately the same ATv,iine instead 
does not significantly alter the results (see Figure lAlY Addition- 
ally, we repeated our analysis replacing the pair-match on metallic- 
ity with a pair-match on stellar mass or stellar surface mass den- 
sity. These latter quantities are global measurements from the total 
SDSS photometry, whereas metallicity is measured from the cen- 
tral 3" of the galaxies. Once again, our results were generally found 
not to be sensitive to exactly which properties we used in the pair- 
matching, presumably by virtue of the observed strong relations be- 
tween these parameters. Where differences were found they were 
consistent with our expectation that dustier galaxies have higher 
metallicities and therefore redder intrinsic SEDs, therefore a fail- 
ure to match on metallicity results in a slight bias in the attenuation 
curves as discussed in the previous subsection. 

We note that no correction is applied for the shifting of the 
photometric bands with redshift (K-correction). We attempted to 
correct each SED ratio for this shifting, using an iterative approach 
in which the previously derived attenuation curve was used to de- 
termine the small corrections to apply to each flux ratio. However, 
this was found to make no difference to the results. The smoothness 
of the attenuation curves compared to individual galaxy SEDs, the 
similar redshift of both members of the pairs, and the relatively 
small redshift range covered by the majority of the galaxies (70% 
with 0.03 < 2 < 0.09) likely all help to minimise the effects. 
For the final derivation of the attenuation curve in Section [6] the 
K-correction is irrelevant as all pairs are combined at their correct 
rest-frame wavelength. 



3 THE DATASETS 

In order to measure the dust attenuation curves from the UV to the 
NIR, in a wide range of star-forming galaxies, we combine three 
of the largest modern day astronomical datasets: the Sloan Digital 
Sky Survey (SDSS) in the optical; the UKIRT Infrared Deep Sky 
Survey (UKIDSS) in the near-infrared (NIR); the Galaxy Evolution 
Explorer (GALEX) in the ultra-violet (UV). 

All fluxes are corrected for Galactic extinction using the val- 
ues provided in the r espective catalogs which are based upon the 
Schl egel et all Il 998) d ust e mission maps and the MW extinction 
curve of ICardelli et al. ( 1989), w ith the sma ll correction in the op- 
tical wavelength regime of lO'Donnelj dl994l) . The GALEX catalog 
provides E(B-V) which we convert into extinction for the UV fluxes 
using the MW extinction curve. 

FigureQ]shows the distributions of redshift and physical prop- 
erties for the galaxies used in this paper. 



3.1 The SDSS spectroscopic survey 

The SDSS is an optical photometric and s pectroscopic survey of 
local galaxies. The final data release (DR7. lAbazaiian et al. 2009) 
includes 9380 sq. degrees with spectroscopic coverage, targeting 
nearly 10 6 galaxies with Petrosian r-band magnitudes < 17.77. 
The spectra have good signal-to-noise and moderate resolution 



which allows the deconvolution of stellar continuum from neb- 
ular emission that is crucial to accurately measure the Ha and 
H/3 line strengths. Such a deco nvolution has been ca rried out 
by iBrinchmann et alj d2004l) and iTremonti et~al d2004l) and the 
emission-line measurements are available onlindj. The SDSS pro- 
vides the u, g,r,i,z aperture photometry used in this paper, with 
central effective wavelengt hs (A e g) of 3546 , 4670, 6156, 7471 an d 
8918A respectively dSchneider et al.||l983l ; iFukugita et ai]|l996l) . 
The SQL query used to extract the parameters from the SDSS Cat- 
alog Archive Server is provided in Appendix |B] A total of 797026 
unique galaxies are extracted from the catalogue, with 552577 that 
have spectral per-pixel signal-to-noise ratio SNR> 6 in the g-band, 
are in the SDSS-MPA database and have "best" magnitudes within 
the formal survey requirements (14.5 < m r < 17.77), where m r 
is the r-band Petrosian magnitude corrected for Galactic extinction. 

Because we are interested in the ISM properties of star- 
forming ga laxies, we use the [N Il]/Ha and [O IIl]/H/3 emission 
line ratios (Bald win et al.lll98ll BPT ratios) to remove galaxies 
which show signs of an Active Galactic Nucleus (AGN), select- 
ing only those galaxie s that lie below the demarcation line of 
iKauffmann et alj I2003T) . We also select only those galaxies with 
all 4 lines measured at >3<r confidence, after scaling the formal er- 
rors by the scaling factors suggested on the SDSS-MPA webpages. 
This results in 140737 galaxies. 

The parameters that we extract from the SDSS catalogues are: 

Emission line dust (tv.iine): the Ha and H/3 line fluxes, mea- 
sured after accurate subtraction of the underlying stellar contin- 
uum, give the total dust attenuation suffered by the emission lines 
at 5500 A. We assume th e dust attenuation curve for emission 
lines of Iwild e TalJ d201ll) : Ty.iine = 3.226 ln(/ H a// H ^)oba - 
3.226 ln(/Ha//H^)int- This was measured through a comparison 
of the luminosity of mid-IR and optical emission lines. In the case 
that the dust free intrinsic line ratio ((fu a / fup)int) is the same 
for both galaxies in a pair, this term cancels when the ratio of two 
galaxy SEDs is calculate^ 

Stellar mass (M*): measured from the 5-band SDSS Petrosian 
photometry in units of solar mass using a Bayesian fit to a li- 
br ary of stochas t ic sta r formation histories similar to those used 
in lGallazzi et alj ( f2005t) are publicly availably Stellar masses are 
calculated from both total and fibre photometry. 

Star formation rate: calculated from the dust corre c ted H a lu- 
minosities, using the conversion factor of Ken nicuttl d 19981) and 
Tyjine given above. We do not use the SFR of Brinchman n et al.1 
(2004) as this assumes a different dust attenuation curve. The pre- 
cise conversions assumed to calculate SFR from Ha luminosity do 
not affect our final results. Specific star formation rate (tps) is cal- 
culated as SFR/M*flb- It was found that ips calculated from the 
5-band SDSS photometry alone was not a good enough estimate 
of the intrinsic spectral shape of the galaxies (see Appendix lAl and 
Figure lA2t . 



° http://www.mpa-garching.mpg.de/SDSS 

6 The intrinsic line ratio depends upon T e ,n e and Z and therefore varies 
slightly with galaxy type and dust content. Over most of parameter space 
relevant to galaxies, this causes a variation in true optical depth for a given 
observed line ratio at a level of less than 10%. However, for massive, high 
metallicity galaxies with low dust contents, optical depths can be up to a fac- 
tor of 2 lower than predicted from the observed line ratio assuming a stan- 
dard intrinsic ratio of 2.87. This effect can not bias our attenuation curves, 
but may act as an additional source of noise in certain bins. 

7 http://www.mpa-garching.mpg.de/SDSS/DR7/Data/stellarmass.html 
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Figure 1. Distribution of basic galaxy properties for the full sample (full line) and UV sample (dotted line). Individual panels show (a) redshift; (b) total stellar 
mass measured from SED fitting to the 5-band SDSS photometry; (c) gas-phase metallicity in the central 3"; (d) stellar surface mass density, calculated from 
the physical size measured in the z-band and total stellar mass; (e) specific star formation rate (ips) measured within the central 3"; (f) ratio of Hce to H/3 
luminosity; (g) axis ratio, measured in the r-band; (h) r-band Petrosian radius in arcsec. 



Gas phase metallicity: measured fro m a combinat i on o f the 
strongest optical emission line fluxes bv lTremonti et al.ld2004l) . us- 
ing the method of Chariot & Longhetti (2001). 

Stellar surface mass density (p*): calculated as M*/(27rr 2 ) 
where r is the physical size in kpc of the radius which contains 
50% of the z-band Petrosian flux. The z-band probes the older stel- 
lar population and thus gives greater weight to galaxy bulges than 
a bluer band would do. 

Size: we adopt the r-band Petrosian radius as a measure of the 
size of the galaxies. This is the radius at which the ratio of local 
surface brightness to mean surface brightness is 0.2. The r-band 
probes both young and old populations with almost equal weight in 
star forming galaxies, and therefore provides a good average mea- 
sure of galaxy size. 

Inclination (b/a): the exponential profile axis ratio, measured in 
the r-band which has the highest SNR of all red SDSS bands. The 
majority of our star-forming sample are expected to be disk galax- 
ies, and therefore the exponential profile fit is the most appropriate. 
A detailed study of the link betw een b/a and t he p hysical shape of 
galaxies in the SDSS is given bv lMaller et all {2009). 

Recent star formation history: we remove all gal axies with 

ost-s tarburst spectral features using the method of IWild et all 

2007t) . These galaxies have unusual SED shapes which may cor- 
relate with dust, thus causing problems during our pair-matching 
procedure. 

A feature of SDSS is that the spectra are observed through 3" 
diameter optical fibres, therefore probing a small central fraction 
of the total light of the galaxies. For the main focus of this work, 
this "aperture bias" is irrelevant as we use the spectroscopic infor- 
mation primarily to identify relative dust content. Where necessary, 
i.e. when we compare fibre and photometric quantities, we use pho- 
tometry extracted in ~ 3" apertures. 



3.2 The UKIDSS Large Area Survey 

We match the SDSS catalogue resulting from the previ- 
ous subsection with the seventh dat a release of the ongoing 
UKIDSS Large Area Survey (LAS, lLawrence etail 120071) to 
obtain NIR Y, J, H, K photometry for the galaxies (A c ff = 
1.0305, 1.2483, 1.6313, 2.2010 Mm). UKIDSS u ses the UKIRT 
Wide Field Camera (WFC A M. ICasali et alj|2007l) . the photomet- 
ric sy_s^ein_J^_^SCTibed in Irlewett et al.l d2006l) , t he calibration 
Hodgkin et ail (2009) and the science archive in lllamblv et al.l 



(2008). We locate all objects with SDSS-UKIDSS cross-matches 
within 2" and select the nearest neighbour if more than one match 
exists. We exclude all galaxies for which photometry has been "de- 
blended". Problems with the UKIDS S pipeline m ean that the flux 
in such objects is overestimated dHill et alj[2O10h . The SQL query 
used to extract the parameters from the UKIDSS online database, 
and cross-match the SDSS and UKIDSS catalogues, is provided in 
Appendix [B] The UKIDSS-SDSS cross match results in our pri- 
mary optical-NIR sample of 22902 star-forming galaxies. 



3.3 The GALEX All-Sky Imaging Survey 

There is considerable interest in the shape of the attenuation curve 
in the rest-frame UV wavelength range, particularly for high- 
redshift galaxy studies. Therefore, we build a second, smaller 
galaxy sample in which we match our optical-NIR gal axy sample to 



the G ALEX All-Sky Imaging Survey catalogue (AIS, Martin et al 



2005). The SQL query used to query the GALEX-MAST databasi 
is shown in Appendix 151 We use the table xSDSSDR7 to identify 



http://galex.stsci.edu/casjobs/ 
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the closest (rank 1) galaxies within 2" jBudavari et a l. 20090, and 
retain only those galaxies with both NUV (A c ff —2267 A) and FUV 
(A e ff =1516A) fluxes measured at greater than 3 a confidence. This 
second UKIDSS-SDSS-GALEX sample contains 15305 galaxies. 



3.4 Aperture photometry 

The precision to which extremely large modern surveys allow us to 
measure stacked quantities demands the use of extreme care in the 
handling of the data if tiny systematics are not to dominate the final 
result. While our pair-matching methodology eliminates absolute 
offsets in photometric scales between different bands and datasets, 
relative offsets can cause noticeable errors. For example, radial gra- 
dients in the stellar populations and dust contents of galaxies mean 
that the use of different photometric apertures at different wave- 
bands would result in incorrect results. Fortunately both the SDSS 
and UKIDSS catalogues provide sets of aperture photometry mea- 
sured with various circular aperture sizes. For each galaxy we inter- 
polate these annular flux densities onto a grid fixed relative to the 
galaxy's r-band Petrosian radius (i? pc t). This allows us to calcu- 
late dust attenuation curves at fixed cumulative fractions of R pe t, 
namely 25%, 35%, 50%, 70%, 90% and 100%. The large PSF of 
GALEX (~5" FWHM) prevents a similar such analysis of the UV 
data, for which we are restricted to a single "total" magnitude. 

Annular aperture photometry is provided in the SDSS cata- 
log at radii of 0.23, 0.68, 1.03, 1.76, 3, 4.63, 7.43 and 11.42 arc- 
sec, which we convert into cumulative fluxes. Similarly, cumulative 
aperture photometry is provided in the UKIDSS catalog at radii of 
0.5, 0.71, 1, 1.41, 2, 2.83, 4, 5, 6, 7, 8 and 10 arcsec. For each 
galaxy, we linearly interpolate both sets of cumulative fluxes onto 
the grid of fixed fractions of R pc t given above. Although linear 
interpolation between measured apertures is not sufficiently accu- 
rate to obtain absolute magnitudes, it is entirely sufficient for our 
purpose of obtaining accurate dust curves. 

The use of fixed-aperture photometry does not account for a 
changing point-spread-function (PSF) as a function of band or sur- 
vey. However, the atmospheric seeing conditions under which the 
SDSS and UKIDSS photometry were obtained were generally ex- 
cellent, and the PSF is typically smaller than 0.35J? pc t (the small- 
est aperture size that we analyse in Section [5} which is >1.5" for 
galaxies in our sample. The small number of very extended galax- 
ies with their Petrosian radius larger than the maximum apertures 
extracted from the catalogs were removed from the sample (see 
Section |3~T| for a definition of Petrosian radius). 



4 THE ATTENUATION CURVES 

To begin with we focus on measurements of the optical-NIR at- 
tenuation curves, for which the data is of a uniform and excellent 
quality, and the data sample larger. In Section IB31 we study the re- 
duced UV-NIR sample of galaxies. 

We split the galaxies into two independent samples with low 
and high stellar surface mass density (/i*), at the position of the 
bimodality obs erved in the loca l galaxy population (i* — 3 x 
10 8 M Q kpc" 2 fcauffmann et alj 120031) . We then define several 
sub-samples in bins of log 10 ^s and b/a. Table Q] presents the de- 
tails of the samples. The measured attenuation curves are given in 



9 http://galex.stsci.edu/doc/CASJobsXTutorial.htm 
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full for each sample and aperture in an online table, Table|2]shows 
an extract. 

The cut in /x* effectively separates bulge-less galaxies from 
galaxies with (even small) bulges. Comparing o ur \i* cut with re- 
sults from the bulge-to-disk decomposition of iGadottil (120090 we 
find that while some pseudo-bulges have 1 x 10 s < fj,* < 3 x 10 8 
and therefore fall into our bulge-less sample, classical bulges and 
ellipticals all have (j* > 3x 10 8 . Radial gradients may be expected 
to be more prevalent in bulge+disk systems, and the presence of a 
bulge may significantly affect the attenuation of lig ht from galaxies, 
compared to pure disk systems dPierini et al.ll2004 . We therefore 
choose to treat these two samples independently. 

Figures|2]and|3]show the attenuation curves for each bin. Dif- 
ferent colour curves indicate different aperture sizes. To facilitate 
comparison between panels, the black line shows the dust curve 
derived from the whole sample, and is the same in all panels. Some 
qualitative comparison between the different attenuation curves is 
instructive. Firstly, we note that in most panels the overall normal- 
isation of the curves drops with decreasing aperture size, which is 
due to radial gradients in dust content. The inner radii are dustier, 
and therefore a larger At is measured. We see that radial gradients 
are more noticeable in galaxies with higher ips. We will present a 
quantitative study of the radial gradients in Section l5~3l 

Secondly, the overall normalisation of the curves in each panel 
depends on the range in dust content of the galaxies in the input 
sample, simply because of the method used to build the curves. A 
larger overall Axv, C ont means a greater range of dust contents in 
the input sample. Both samples show a steady decrease in ATv lCon t 
with increasing b/a and increase in ATv,cont with increasing ips- 
This implies that high ips and more inclined galaxies exhibit a 
larger range of ATv,cont- This pure selection effect does not af- 
fect our results, which depend only on the shape of the attenuation 
curve, or its amplitude relative to another dust indicator. 

It is noticeable from Table Q] that the bins of ips and axis ra- 
tio are not entirely independent: when binning by ips, bins with 
lower ips have lower median axis ratios, and similarly when bin- 
ning by axis ratio. Indeed, when investigating the whole dataset 
there is a small but significant positive trend between ips and axis 
ratio, which cannot be a "fundamental" property of galaxies but, as 
we shall show later in the paper, is likely an effect of dust. A bias 
in an observational property such as this one can cause an intrin- 
sic property to correlate with dust content, and ultimately bias the 
shape of the attenuation curve. In this case it shows the importance 
of pair-matching galaxies in b/a. 

In the following section, we begin our quantitative study of 
the attenuation curves of the binned samples by parametrising them 
using a simple analytic form, which allows for easy visualisation of 
changes of shape and amplitude with galaxy properties. 



5 TRENDS OF ATTENUATION WITH GALAXY 
PROPERTIES 

As discussed in Section Q] the shape of the dust attenuation curve, 
and strength relative to the ratio of Ha to H/3 luminosity, can tell us 
about the dust properties and dust-star geometry within galaxies. In 
order to visualise trends in the attenuation curves presented in Fig- 
ures|2]and|3]we must first parametrise the curves. After some exper- 
imentation, a broken power-law was found to provide an adequate 
representation of the optical-NIR data, with easily interpretable fit- 
ted parameters. A single power-law was unable to fit the whole 
wavelength range, due to a change in the slope between the optical 
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Table 1. Number of pairs, median redshift, median ip„$> median b/a, and mean Atvjine for each SDSS+UKIDSS sample. 
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AT V,line 
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All galaxies 
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0. 
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0. 


.61 
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< log^s/yi^ 
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0.072 
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Table 2. Averaged attenuation curves, normalised to zero at K (7a, Eqn.|7J, for each sample studied in this paper [full table available online]. 



Sample 


Aperture 


u 


9 


r 


i 


z 


Y 


J 


H 


All galaxies 


0.25i?. pct 


-0.493 


-0.362 


-0.261 


-0.209 


-0.161 


-0.123 


-0.074 


-0.042 




0.35Rpet 


-0.473 


-0.349 


-0.250 


-0.200 


-0.153 


-0.116 


-0.068 


-0.038 




0.50_R pct 


-0.455 


-0.336 


-0.239 


-0.192 


-0.146 


-0.110 


-0.062 


-0.035 




0.70.Rp C t 


-0.439 


-0.325 


-0.231 


-0.186 


-0.141 


-0.105 


-0.056 


-0.033 




0.90i? pct 


-0.431 


-0.319 


-0.227 


-0.183 


-0.139 


-0.101 


-0.054 


-0.032 




1.00i?. pct 


-0.428 


-0.316 


-0.225 


-0.182 


-0.138 


-0.100 


-0.054 


-0.032 




3" 


-0.484 


-0.355 


-0.253 


-0.202 


-0.154 


-0.119 


-0.069 


-0.038 



and NIR wavelength regimes. We fit this function using a least- 
squares minimisation routine, weighting each data point equally: 

Tx = -At a I— j +At k A<A c 

= -At Ac (^-pj +At k A>A c (9) 

where the free parameters describe the wavelength at which the 
shape of the attenuation curve changes between the optical and the 
NIR (A c ), the mean amplitude of the attenuation curve at this wave- 
length (Ata c ), the mean amplitude of the attenuation curve at the 
longest wavelength point (Atk), the power-law slope of the attenu- 
ation curve in the optical (s op t) and in the NIR (s nn ). In this section 
our aim is only to measure the slope of the curves for the purpose 
of visualisation. Therefore, we fit only the u, g, r, i, Y, J, H and K 
bands, excluding the z-band because of its proximity to the break 
in slopes at A c . For the same reason we do not constrain the am- 
plitudes of the two powerlaws (Ata c ) to be equal. In practice the 
amplitudes of the two functions are similar. Note that by construc- 



tion the form of this function is the same as that derived in Eqn.[8] 
with Q\ (x X s . 

Bootstrap errors on each of these parameters are estimated by 
recalculating the geometric mean, using random resampling with 
replacement of the galaxy-pairs that contribute. We randomly se- 
lect 100 samples of the same size as the original dataset, and repeat 
the least squares fitting on each sample. The 16th and 84th per- 
centiles on each parameter are taken to be the lcr errors. We veri- 
fied that increasing the number of samples to 500 does not change 
the resulting estimated errors. 

Figures [2] and [3] show that even in samples exhibiting strong 
radial gradients, the attenuation curves are close to converged be- 
tween 90% and 100% of R pc t- The fitted parameters at 90 and 
100% of i?pct are equal within the errors, but the errors at 90% 
are smaller and trends tighter. This is possibly caused by increas- 
ing sky subtraction errors at large radii in the UKIDSS photometry. 
Therefore, unless otherwise stated, we measure the galaxy-wide at- 
tenuation curves at 90% of i? pe t. 

The fit is performed at the median rest-frame of the galaxies in 
the sample (the median redshift is given in Table [T). As mentioned 
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Figure 2. Averaged dust attenuation curves, normalised to zero at K (Eqn.|7J. Each panel shows galaxies in bins of tpg, top: with log /i* < 8.3, and bottom: 
with log fi* > 8.3. Note that the overall normalisation of the curves in each panel depends on the range in dust content of the galaxies that make up each 
sample, and is thus somewhat arbitrary (see text of Section |4j . Different colour lines indicate different physical aperture sizes, from red to purple: 25, 35, 50, 
70, 90 and 100% of the Petrosian r-band radius. In some panels the curves show a progressively lower normalisation with decreasing aperture size. This is 
caused by larger dust contents in the central regions of the galaxies. The dotted black line shows the attenuation curve when a 1.76" radius angular aperture is 
used, close to the aperture through which the SDSS spectra are taken. To aid comparison between panels, the full black line is the attenuation curve measured 
at 90%i?pct using all galaxies, and is the same in each panel. 



in Section |231 iteratively K-correcting the individual flux ratios be- 
fore combining them into the dust attenuation curve makes no sig- 
nificant difference to the results and was therefore discarded from 
the final analysis. 

The resulting fitted parameters and errors are given in Table [3] 
together with comparison numbers for curves in the literature mea- 
sured in exactly the same way as our data. As expected, the fitted 
Atk values are close to zero, with their exact values imposed by the 
constraint that the NIR slope is a perfect power-law. In the follow- 
ing subsections we present the trends of attenuation curve shape, 
radial gradients and difference between nebular line and stellar con- 
tinuum attenuation as a function of ifrs an d axis ratio, for both high 
and low /i* galaxies. 

5.1 The slope of the attenuation curve in the optical 

Figure |4]presents the optical slopes of the attenuation curves (s op t) 
as a function of tps and b/a. The dotted horizontal lines indicate 
the slopes of the MW, COO, LMC and SMC-bar curves measured 



using the same method. There are significant trends in optical slope 
with galaxy properties. In both high and low /i* samples the curve 
flattens slowly with decreasing axis ratio, meaning greyer curves 
in edge-on galaxies and steeper curves in face-on galaxies. This is 
the first clear example of global geometry affecting the apparent 
attenuation curves of galaxies. 

The high n* sample exhibits steeper optical slopes than the 
low [i* sample, and significantly steeper even than the Milky Way 
extinction curve. This is intriguing as attenuation curves resulting 
from mixed stars and dust in galaxies are generally expected to be 
flatter than the underlying extinction curves. Steeper optical extinc- 
tion curves do exist locally in the SMC and actively star-forming re- 
gions in the MW, possibly caused by the modification of grain com- 
position. However, it would be surprising if massive, metal rich, 
bulge dominated galaxies have grain compositions closer to active 
dwarf galaxies than to the MW. Steeper apparent slopes could be 
caused by a mismatch in the stellar population of the galaxy pairs. 
However, this requires the more-dusty galaxies to have an intrin- 
sically older and therefore redder stellar population than the less- 
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Figure 3. Averaged dust attenuation curves, normalised to zero at K (see equation |7J. Each panel shows galaxies in bins of axis ratio (a/b), top: with 
log fi* < 8.3, and bottom: with log fi* > 8.3. Lines as in Figure|2] 
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Figure 4. Power-law slope of attenuation curve in the optical wavelength range, as a function of rpg (left) and axis ratio (right), for low fi* (open circles) 
and high /i* (tilled circles) galaxies. The horizontal lines indicate the equivalent slopes for the MW (dotted), COO (gray, full), LMC (dashed) and SMC-bar 
(dot-dash) curves. 



dusty galaxies, despite having the same observed ips, metallicity 
etc., which is opposite to the slight trend with inclination that we 
observe in the data. 

It is possible that the steeper slopes are a result of a differen- 
tial effect related to the distribution of dust around stars of different 
ages. In ordinary star-forming galaxies, younger (older) stars con- 
tribute a larger fraction of the blue (red) optical light. This is dif- 
ferent from low-mass starburst or old elliptical galaxies, in which 
the majority of light at all wavelengths comes from young or old 
stars. If young, blue stars suffer greater dust extinction than old, 
red stars, then we will observe an apparently steeper attenuation 
curve in star-forming galaxies than the underlying extinction curve. 
This effect is the common explanation for the difference in attenu- 
ation observed in emission lines and galaxy continua (see Section 
15.4b . We observe that the slope is steeper in galaxies with signifi- 



cant bulges. This could occur if the older stellar population in the 
bulge is less obscured by dust than the younger stellar population in 
the disk. Slopes closer to the extinction curves observed in edge-on 
galaxies could arise as the disk obscures a greater fraction of the 
old stellar population in the bulge. 

If the steeper slopes are a result of the relative geometry of 
stars and dust, we would expect a weaker effect in the smaller aper- 
tures where the light comes mainly from older stars in the bulge, 
especially in face-on galaxies: at 0.35i? pc t we find s op t = 1.07 for 
fe/a=0.35, and 1.31 for b/a=0.88. This is 3% and 22% smaller 
than the values for s op t at 0.9R pct , qualitatively in agreement with 
the hypothesis that the steep curves arise due to a differential ef- 
fect whereby different stellar populations are attenuated by differ- 
ent dust columns. 

One final possibility is presented by radiative transfer models 
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Table 3. Parameters of broken power-law fits to attenuation curves for each SDSS+UKIDSS sample (see Eqn.[9|l. Measurements of the shape and amplitude 
of the attenuation curves are made from circular aperture photometry measured at 90% of the r-band Petrosian radius. The line to continuum ratio is measured 
using a 1.76" radius fixed angular aperture, to best match the emission line measurement through the SDSS fibre. Note that the slopes and t ratios are 
fundamental quantities, whereas the At values are sample dependent. 
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a: E(B -V) sta . IS = 0.44£fB-Vj gas (Calzetti 1997bL l200lh converted into a ratio of optical depths assuming the emission line attenuation curve from lWild etafl 

Honl). 
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of dust in galaxies. These can also exhibit steeper optical attenu- 
ation curves than the input extinction curves, in the optically thin 
limit and with cer tain geometric configurations dTuffs et alJl200l 
IPierini et al]|2004l) . even before effects for different stellar popula- 
tions are included. Ultimately, the cause for the steep optical slopes 
must be solved by comparison with such models. This is beyond the 
scope of the present paper (see Chevallard et al., in preparation). 



5.2 The slope of the attenuation curve in the NIR 

In Figure|5]we show the slope of the attenuation curve in the NIR, 
as a function of ips an d axis ratio. Within the errors we find no 
significant overall deviation from the MW extinction curve with a 
slope of 1.6, and no clear trends with galaxy properties are evident. 
For the LMC and COO curves we find NIR slopes of 1.7, almost 
identic al to the MW curve. The SMC-bar curve of iGordon et al.l 
■2003h exhibits a strong bump in the NIR that prevents us from fit- 
ti ng a robust powerlaw slope. Inspection of the individual curves 
in IGordon et al] ,20031) reveals that a significant deviation from a 
power law slope is only evident in one out of the four stars that con- 
tribute to the average curve. We therefore conclude that this bump 



may result from photometric errors and small number statistics, and 
we therefore do not compare our results with the SMC-bar curve in 
the NIR wavelength regime. 



The constancy of the NIR slope of extinction laws along mul- 
tiple lines of sight to stars in the MW is an important feature of 
the MW extinction curve, implying a constant size distribution of 
the l argest particl es throughout both diffuse and dense clouds of 
dust (Mathis 1990). The same value found here in external galaxies 
leads to the conclusion that the size distribution of the largest dust 
particles is universal. It additionally implies that the factors affect- 
ing the attenuation of light in galaxies (dust-star and global geome- 
try, scattering into the li ne-of- sight) are unim portant at long wave- 
lengths. Recent work by Stead & Hoard d2009l) suggests a slightly 
steeper MW NIR slope of 2.14 by comparison of UKIDSS data 
with the Besancon Galactic model. If this is the case, this would in- 
stead imply some degree of flattening in the NIR attenuation curve 
of external galaxies. It is clear that further work on measuring the 
shape of extinction curves in the MW and local galaxies would be 
helpful to understand the dust distribution in more distant galaxies. 
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Figure 5. Power-law slope of attenuation curve in the NIR wavelength range, as a function of tp,$ (left) and observed axis ratio (right), for low (j,* (open 
circles) and high (i* (filled circles) galaxies. For clarity the high fi* symbols have been offset slightly to the right. The dotted line indicates the equivalent 
slope for the MW extinction curve. 
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Figure 6. Radial gradients in continuum optical depth. Relative optical depth at 5500A towards stars contained within 35% (yellow, upwards triangle) and 
50% (green, downwards triangle) of the Petrosian radius of the galaxies, relative to the optical depth measured at 90% of the Petrosian radius, as a function of 
tpS (left) an d a "i s ra ti° (right). The top (bottom) panel shows the results for low- (high-) fi* galaxies. 



5.3 Radial Gradients in total dust attenuation 

In this section we provide a more quantitative analysis of dust ra- 
dial gradients than can be achieved by visual inspection of Fig- 
ures [2] and [3] Because our methodology only provides a rela- 
tive difference in dust content between the less- and more-dusty 
galaxies, we normalise the measured Atv values at all radii by 
the value at 0.9i? pc t. This then allows a direct comparison be- 
tween samples. We measure Atv from the independent power- 
law fits to the dust curves at each different radius. It can be 
shown that, in the case that Tv(r)/Tv(0.9i? pc t) is the same on 
average for both the more and less-dusty galaxies that make up 
the pairs in a particular sample, i.e. within a particular bin the 
strength of the radial gradients is independent of overall dust con- 
tent, then AT V (r)/ATv(0.97?pet)=Tv(r)/Tv(0.9i?p Ct ). However, 
whether this is the case or not does not affect the identification of 
trends in which we are primarily interested. 

Figure [6] shows Atv measured within 0.35 and 0.5i? pe t, rel- 



ative to Atv measured at 0.9i? pc t. In this section only, we fix the 
NIR slope to 1.6 during the fits, because for small apertures the er- 
rors increase substantially on the measured NIR slope. In the left 
hand panels we see a strong trend in the strength of radial gradients 
with tps'- gradients are significantly larger in high ips galaxies than 
low tps galaxies. In the low fi* sample, the trend appears to flatten 
above ips^ —9.5. In the right hand panels, we see no strong trends 
in the strength of the radial gradients with axis ratio, although the 
slightly weaker radial gradients in very inclined galaxies is likely 
a real feature caused by large optical depths. It is noticeable that 
both the high and low fi* samples show the same strength of ra- 
dial gradient. Averaging over all bins in b/a, we find that the dust 
optical depth within 0.35 and 0.57? pc t is 1.16 and 1.09 times that 
within 0.9i? po t, for both samples. This indicates that the presence 
of a bulge does not affect the strength of radial gradients in the dust 
distribution. 

Measuring the trends of radial gradient strength with ips 
is important for understanding the production and destruction of 
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dust relative to ongoing or past star formation. The dust gra- 
dients reported here, and elsewhere in the literature for spa- 
tially resolved studies of small numbers of local galaxies (e.g. 
Muno z-Mateos et al] |2009). are presumably related to metallicity 
gradients in spiral disks. However, because our ips measurements 
relate strictly to the central 3" of the galaxies, and the scale of star 
formation episod es is known to depend on the rate of star formation 
jKennicutH l998). spatially resolved information on star formation 
rates is required before firm conclusions can be drawn. Suitable 
spatially resolved spectral studies of metallicity, dust and star for- 
mation gradients in large numbers of disk galaxies have yet to be 
undertaken. 



5.4 Dust attenuation of nebular emission lines vs. stellar 
continua 

All of the galaxies in our sample have both Ha and H/3 line mea- 
surements which, until now, we have only used to derive relative 
dust content. Here we make use of the line fluxes to compare the 
dust attenuation suffered by the nebular emission lines and stellar 
continuum. The SDSS spectra are taken using a 3" diameter fibre, 
meaning that smaller galaxies, and those at higher redshift, have a 
larger fraction of their total light included in the spectrum. In our 
sample, the median projected fibre radius is around 0.35 -R pe t, but 
with a wide distribution. In light of the radial gradients measured in 
the previous section, this fixed angular aperture restricts the extent 
to which we can compare dust attenuation of lines and continua in 
our sample. To provide the best comparison, in this section only, 
we measure the attenuation curves using a fixed 1.76"radius angu- 
lar aperture (the curves are shown as dotted lines in Figures [2] and 
[3}. This is slightly larger than the 1.5"radius aperture of the fibre, 
however the spectra are taken in worse seeing conditions than the 
imaging, and flux measurements at 1.76"are provided in the SDSS 
catalog thus preventing any small errors introduced by interpolation 
of the radial flux curves. 

Using the same formalism as for dust attenuation of the stellar 
continuum (Section r2.lt . the amount of dust attenuation of nebular 
emission lines in more-dusty relative to less-dusty galaxies is given 
by: 

AXA.Iine = ATv,lineQA,line (10) 

where QA.iine describes the shape of the attenuation curve applica- 
ble for emission lines. Wil d et all d201 ll) compared mid-IR and op- 
tical emission line str engths to show that a two component model 
(Chariot & Fall 2000), with diffuse dust accounting for a fraction 
(fj.) of the optical depth at 5500A, and denser birth-cloud dust ac- 
counting for the remainder, is consistent with the emission line ra- 
tios observed in a wide range of galaxy types: 

QA,!ine = (1 - ft) (A/Ay) - ™ + /t(A/Av) _ ™ (11) 

In IWild et al.1 ( 1201 ll) we measured fj, = 0.4, for the case that 
m — 1.3 and n = 0.7. Ultimately the slopes n and m may also be 
measured through comparison between the attenuation curves mea- 
sured in this paper and models. However, for the purposes of this 
paper it is sufficient to know that an accurate Tyj j ne is obtained from 
the Balmer emission line ratios, as shown by IWild et all d201 lh . 
Eqn.[TT]leads to the dust optical depth for emission lines, as a func- 
tion of wavelength and ratio of Ha to H/3 luminosity: 




where Q H ^,iine - QHa.line = 0.310. 

In Figure|7]we show the ratio of dust optical depth in the emis- 
sion lines to optical depth in the stellar continuum as a function of 
ips and axis ratio. In both samples we find a strong trend with ips 
(left panel), with both low and high /i* galaxies having a line-to- 
continuum attenuation ratio 50% higher at lower ips than at high 
ips- There is also a strong trend with axis ratio, with high axis ratio 
galaxies having a higher line-to-continuum attenuation ratio than 
low axis ratio galaxies. Similar trends with axis ratio are observed 
in both low and high /i* galaxies, but low /i* galaxies are generally 
offset to higher overall ratios. 

These trends can be interpreted in terms of the two different 
dust components in galaxies (diffuse vs. birth-cloud dust as de- 
scribed above), with young and old stars suffering different atten- 
uation. The difference between the high and low n* galaxies may 
relate to the lower overall metallicity of the diffuse ISM in the lat- 
ter sample, leading to lower diffuse dust contents. The slightly non- 
linear trend with b/a for the low n* galaxies follows qualitatively 
that expected for a dust disk of finite thickness, with the optical 
depth of the diffuse ISM varying with the ratio of scale-height to 
scale-length. One further point of note is the equal ratio found at 
high ips and low b/a. Such a coincidence of ratios is not expected 
a priori, as the former property is intrinsic to the galaxies and the 
latter is observer dependent. This may indicate the true limit of op- 
tical observations, where significant numbers of sight lines are lost 
due to high column densities of diffuse dust, either due to long path 
lengths or high dust contents. Such a "skin" effect would lead to 
only the stars closest to the observer contributing to the observed 
integrated light. We will return to all these points in Section|7] 

As in the case of the radial gradients presented in the previous 
subsection, in the case that Tv,iineAtV,cont ls independent of overall 
dust content, it can be shown that ATv,iine/ATv. CO nt=T : v.iine/Tv,cont- 
In Appendix [A] Figure [ATI shows that this is indeed the case, and 
therefore we can write (dashed lines in Figure[7J: 

= 0.75(±0.2) +3.8(±0.6)6/a (law/**) (13) 

TV.cont 

= 1.6(±0.2) + 1.4(±0.4)6/a (high /t*) (14) 

and 

= -6.4(±2.6) - 0.9(±0.3)V> * (low/**) (15) 

TV,cont 

= -8.0(±2.6) - 1.1(±0.3)V> * (high/w) (16) 

ICalzettil (Q997b) measured the same quantity for starburst 
galaxies, finding E(B— V) star s = 0A4E(B— Vj\ ln& . For the emis- 
sion line attenuation curve given in Eqn. [TT] this corresponds to 
tv, line/iv, com = 0.44R v ,ii n e /Rv.cont = 2.08. This is shown as a 
dotted line in Fig. [7] which is close to the value measured for high 
ips galaxies in our sample. 



5.5 The slope of the attenuation curve in the Ultraviolet 

In the previous subsections we used the larger SDSS+UKIDSS 
galaxy sample. We now turn to the analysis of the subsample of 
galaxies with UV fluxes, repeating the pair-matching procedure 
from the start. The results of the optical-NIR fits remain the same 
as in the full sample, albeit with larger errors resulting from the 
smaller sample sizes. In this section we define the slopes for FUV- 
NUV and FUV-u as: 

g > = log(-TAi) -log(-7A 2 ) 

log Ai - log A 2 
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Figure 7. Dust optical depth towards the Balmer emission lines, relative to dust optical depth in the continuum at 5500A, as a function of tps (left) and axis 
ratio (right). Unlike previous figures, the optical depth in the continuum has been measured using a fixed angular aperture of radius 1.76", to approximately 
match the emission line measurements which are obtained through an optical fibre with 1.5"radius aperture in generally worse seeing conditions. The dotted 
line indicates the line to continuum optical depth ratio from Calzetti 1 1997b), assuming the same attenuation curve for the emission lines as used in our analysis 
(Eqn. lllt . The dashed lines are fit to the data points and are described in the text. 



4 




2000 3000 4000 5000 

Observed wavelength / urn 

Figure 8. The MW (full) and COO (dashed) attenuation curves in the UV- 
optical wavelength range. Overplotted are the UV filter response functions 
from GALEX and the SDSS u-band. The stars (MW) and crosses (COO) 
indicate the attenuation in each of these bands after convolution with the 
filter functions. 



Note that s' would only be equivalent to the power-law slopes (s) of 
the previous subsections (Eqn. [5), if the attenuation in the K-band 
(Atk) were equal to zero. Due to the smaller dataset used in this 
sample, we quote s' rather than s in order to provide measurements 
that are independent of the optical-NIR attenuation curves. At these 
short wavelengths, s and s' are not significantly different. We test 
this by measuring both s and s' from the literature dust curves, nor- 
malising at the rest-frame wavelength of the K-band at the median 
redshift of our samples. We find s' — s ~ 0.03, well within our 
measurement errors. 

To quantify the influence of the strong bump in the MW ex- 
tinction curve on broadband attenuation measurements, we con- 
volve the comparison dust curves with the filter response functions 
(see also the discussion in Section [2~2l l. Figure[8]shows the MW and 



COO dust curves in the UV, with the FUV, NUV and u-band filters 
overplotted. The stars/crosses indicate the filter-averaged attenua- 
tion of each curve in the three bands: clearly the filter- averaged 
attenuation in the GALEX NUV band for the MW curve is less 
than the actual attenuation at 2175A. This figure also shows how 
the slope between the FUV and NUV bands, at median rest-frame 
wavelengths of 1430A and 2138A, probes the depth of the 2175A 
dust feature, whereas the slope between the FUV and u bands mea- 
sures the shape of the underlying dust continuum. 

Table|4]gives the measured UV slopes for all samples, together 
with comparison values for dust curves in the literature. In the top 
panels of Fig. [9] we find that s FUV _ u is slightly higher on aver- 
age in the high n* samples, than the low n* samples. The two 
samples straddle the value measured for the MW and COO curves 
(Spuv-u— 0.6 in both cases). There is some evidence that the UV 
attenuation curve steepens with increasing b/a in the high fj* sam- 
ple, and there is a possible steepening of slope with tps in the low 
/x* sample. 

In the bottom panels of Fig. [9] we find that s PUV _ NUV gen- 
erally lies between the values measured from the MW and COO 
curves. The MW extinction curve, with a strong 2175A dust-bump, 
has Sp UV _ NUV = 0.06, whereas the COO attenuation curve with no 
dust bump has Sftjv-nuv = 0.58, similar to the global (FUV-it) 
slope of this curve (sftjv-u— 0.62). For the high fi* galaxies we 
find that the bump strength increases strongly with increasing b/a, 
ranging from close to the COO curve in face-on galaxies, to close 
the MW curve in edge-on galaxies. The same trend is not observed 
in the low /i* galaxies. A weaker trend of increasing bump strength 
with decreasing tps may reconci le our results with th e weak bump 
observed in starburst galaxies bv lCalzetti et al.l (12000). 

The 2 175 A feature is widely believed to be a pure absorption 
feature, with its strength in extinction curves varying dependent on 
the che mical compo sition and/or physical environment of the dust 
grains JPrainell2003l) . In attenuation curves, the balance of scatter- 
ing vs. absorption could affect the strength of the feature. The sim- 
ilarity in overall strength between the low and high-/i* galaxies, 
which have different average metallicities, could imply that much 
lower metallicities than are present in our sample are required to 
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Table 4. Slopes of the UV attenuation curves for each 
SDSS+UKIDSS+GALEX sample (see Eqn. [T7J), and comparison lit- 
erature curves. 



s' s' 
a FUV-NUV s FUV-u 



H* < 3 x 10'' 





All galaxies 





no_|_0.04 
zo=t 0.03 





c 4 _i_0.02 
04=t 0.03 


-11.00< 


Mife/V" 1 ) <-9.90 





nq+0.14 
UM=t 0.12 





oc+0.09 
JD=I: 0.08 


-9.90< 


log^s/yr" 1 ) <-9.75 





2fi+ ' 17 





44+0.08 
44=t 0.08 


-9.75 < 


logCVs/yr" 1 ) <-9.62 





47+ 0.l4 





coj.0.10 
b8± 0.07 


-9.62< 


logCVs/yr- 1 ) <-9.46 





oq+0. 09 
ZJ3= 0.08 





4K+0.07 
* o=t 0.06 


-9.46< 


logtVs/yr- 1 ) <-9.25 
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ZD=C 0.09 
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log^s/yr" 1 ) <-7.90 





41+O.H 
41=c 0.11 





fi c:+0.07 




0.0< b/a <0.4 





on+0.06 





ec:_l_0.05 
oo=c 0.05 




0.4< 6/a <0.6 





on+0.06 





c fi _i_0.06 
OD=C 0.04 




0.6< b/a <0.8 





in+0. 12 
iy=t 0.08 





cq._i_0.08 
O,5=l: 0.06 
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oo=c 0.12 





coiO.08 
oc>=l: 0.08 



^t* > 3 X 10 s 


All galaxies 
-11.00< log^g/yr" 1 ) <-9.90 
-9.90< log^s/V" 1 ) <-9-75 
-9.75< logW-s/yr- 1 ) <-9.62 
-9.62< log^s/yr" 1 ) <-9.46 
-9.46< log^s/yr -1 ) <-7.90 

0.0< 6/a <0.4 

0.4< b/a <0.6 

0.6< b/a <0.8 

0.8< 6/a <1.0 


QO_|_0.04 
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o-2i±8:il 

0.30±° : « 
0.39±o : ;? 

n 2S+° 10 
u.zo_r og 

n 4S+O.1O 

o.56±g : °? 


0.74±°;°| 

0.88±g;0? 
n fifi+o.io 

n 70+0. 09 

n 72+0-07 
u -' z=c 0.07 
n K7+0.06 
u -°' =t 0.05 
n 67+0.07 

n fi4+0-06 
U - D4=t 0.06 
n 7C+O.O6 

n a«+o.o8 


MW (CCM+OD) 


0.06 


0.64 


COO 


0.58 


0.62 


LMC 


0.27 


0.71 


SMC-bar c 


1.24 


1.24 



reduce the strength of the feature through chemical composition 
effects. The samples studied here are dominated by luminous and 
relatively massive galaxies. The COO curve is well known for the 
lack of 2175A feature, which may be due to dust grain destruction 
in the extreme physical environment. Again, although there is a hint 
in the data that the strength of the 2175A feature weakens with in- 
creasing ips, the samples studied here are dominated by ordinary 
star-forming galaxies, not extreme starbursts. 

The trend with axis ratio in the high /i* galaxies is interesting, 
as is the fact that the same trend is not seen in the low fj,* galax- 
ies. The lack of overall difference in strength between the high and 
low /_* samples implies that metallicity effects are not important 
here, which leaves the possibility of geometric effects. If dust is 
distributed uniformly with the stars, then by increasing the incli- 
nation of the galaxy, the amount of interstellar dust that star light 
passes through on its way to the observer increases, which will in- 
crease the optical depth of absorption by the grains responsible for 
the 2 175 A feature. However, in the case of high /1* galaxies which 
have significant bulges, axis ratio is not as cleanly linked to incli- 
nation as in the low fi* galaxies. This makes our results difficult to 
interpret without further image decomposition that is not possible 
with GALEX data. 

We conclude that a 2175A dust feature is present in most types 
of luminous galaxies, except possibly at high ips, but its observed 
strength in attenuation curves can be modified by geometric ef- 
fects. How these geometric effects act will require better spatially 
resolved observations and more detailed evaluation of dust models. 



The implications of a significant 2175A dust feature extends 
beyond interest in the composition of dust grains, as it impacts 
significantly on the observed UV spectral slope of galaxies and 
thus estimates of the s t ar formation rate de nsity at high red- 
shift .Meureretailll999l; iBouwens et alj|2009), and the metallic- 
ity and dust contents of high redshift gal axies Meurer et alJI 19971 : 
IBouwens et alj201ol ; lDunlop et al.ll__.01 lb . Predictions for the exact 
strength of the effect are difficult given the wide variety of bands 
used to measure the UV spectral slope (commonly referred to as /3), 
and the lack of constraints on the stellar population of high redshift 
galaxies. However, when the reddest band used to measure f3 lies 
close to 2175A as in lBouwens et al.l d2010h . incorrectly assuming a 
dust curve with no 2175A feature would lead to an underestimation 
in the amount of dust, leading to the incorrect conclusion that the 
stellar populations are extremely metal poor. 



6 GLOBAL FIT OF THE DUST ATTENUATION CURVE 

In the previous section, we used a combination of a parametrised fit 
and direct colour measurements to study the shape of the attenua- 
tion curves for independent galaxy samples binned by ips and axis 
ratio. Parametrised fits of binned samples are important for robustly 
identifying systematic changes, but they do not provide an "attenu- 
ation law" that can be used at all wavelengths on all datasets. In this 
section, we provide an empirically derived curve to describe the at- 
tenuation by dust at arbitrary wavelength, as a function of ips and 
axis ratio, suitable for any observational dataset. Where ips and 
axis ratio are unknown, and no informed guess can be made, the 
"typical" curve may be defined at the central points of ips^ —9.5 
and axis ratio 6/a~ 0.6. 

In this section the data are not binned in ips or b/a, so that 
each galaxy pair can be included at its rest-frame wavelength and a 
smooth function with ips and b/a can be derived. We use a multi- 
ple smoothly broken power-law to describe the shape of the curve 
.Guidorzi etaT]|2009l) . 

In order to derive an average attenuation curve shape (Q\) 
without prior binning of the data, we require an estimate of the op- 
tical depth in the continuum (ATv, con t) for each pair, a quantity not 
known a priori. We therefore use the emission line optical depth, 
taking into account the variation in line to continuum dust opacity 
(ATvjine/AT V ,cont) with ips and b/a observed in Section [54] and 
the radial gradients in overall dust content observed in Section [573l 

With only 2 broad bands in the UV, an accurate description 
of the 2 175 A region of the attenuation curve is not possible. While 
there is a definite change in slope of the curves, indicative of the 
2175A feature observed along most lines-of-sight in the MW, we 
are unable to constrain the precise shape of this feature. We have 
therefore chosen to define a break in the attenuation curve at the 
central wavelength of the feature observed in the MW, but inter- 
polate smoothly over the broad band observations. This will not 
cause significant errors in the correction of broad band fluxes for 
dust attenuation. However, when using our empirical curve to cor- 
rect spectral observations, or narrow band fluxes, readers should be 
aware of the limitations of our observations. 

The curve is composed of four power-law functions with ex- 
ponents S[!_ 4 ] smoothly joined with a smoothness parameter n. We 
define the position of the three break points (Ab[i_3], e ff) at 0.2175, 
0.3 and 0.8/im. The first is motivated by the central position of the 
dust bump in the MW curve. The second break point is positioned 
to produce a smooth transition between the steeper optical and shal- 
lower NUV slopes observed. The final break position is motivated 
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Figure 9. Power-law slope of the attenuation curves in the UV wavelength range, as a function of ips (left) and observed axis ratio (right). In the top panels we 
plot FUV-u band slope. The lower panels show the FUV-NUV slope. For clarity in the lower panels the high /j,* symbols have been offset slightly to the right. 
The horizontal lines mark the corresponding slopes for the MW (dotted), COO (gray, full), LMC (dashed) and SMC-bar (dot-dash) curves after convolving with 
the filter response functions. The FUV-m and FUV-NUV slopes of the SMC-bar extinction curve are 1.24, which is beyond the axis range. 



Table 5. Fitted constants for Eqn.[l8] which defines the shape of the dust attenuation curve as a function of tpg and axis ratio. 



Sample 



/i*<3x 10 1 
At*> 3 x 10 1 



CI 


C2 


C3 


C4 


C5 


C6 


C7 


C8 


C9 


0.15 
0.20 


0.00 
0.90 


0.20 
0.10 


0.70 
1.10 


0.00 
0.30 


0.40 
-0.20 


1.10 
1.30 


0.40 
0.60 


-0.10 
0.10 



by the results of our parametrised fits in Section|5] The sampling of 
the curves is too sparse to fit for the smoothness parameter n, and 
therefore we fix this to be 20 which results in a sufficiently smooth 
curve. 

Following the formalism set out in Section l2~2l the shape of 
the attenuation curve is given by: 



Qx 



1 

~'N 



A 



A 

Abi 



n .so 
+ 



A 

Ab2 



A 

Ab3 



-1/n 



(18) 

where N is the normalisation, defined such that the curve is unity 
at 5500A, Qv = 1; n — 20 defines the smoothness of the breaks. 
Note that due to the smoothing, the exponents of Eqn. Q~8] are not 
all trivially related to the slopes measured in the previous section. 
The effective position of the break points are converted into 



the A b [i-3] values according to: 



Abi 

Ab2 
Ab3 



(Ab2,eff) S 



(A 



b3,eff) 



l/*3 



l/s 



(19) 
(20) 

(21) 



with the effective positions of the breaks as motivated above. 

The power-law slopes are allowed to vary with both axis ratio 
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Figure 10. Dust attenuation curve shape {Q\), as given by Eqn. 1181 as a function of tps (top) and axis ratio (bottom) for galaxies with low fjt* (left) and high 
/i*(right). Attenuation curves at values of log 10 ips= [—9.9, —9.7, —9.5, —9.3, —9.1] and b/a= [0.2, 0.4, 0.6, 0.8, 1.0] are shown as different style lines, 
as indicated in the top right of each panel. When variation with ips (b/a) is shown, b/a (log 10 ips) is held constant at 0.6 (-9.5). 
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Figure 11. The observed shape of the dust attenuation curves in the binned samples (with UV data: crosses; without UV data: stars [i.e. normalised version 
of Figs.|2]and[3)), compared to the empirical dust curve shapes derived from the entire (unbinned) samples (line and diamonds). The median ipg and b/a of 
the binned samples have been used to extract the appropriate attenuation curve from Eqn. 1181 these are given in the bottom left of each panel. The left (right) 
panels show the low (high) fi* samples. The top (bottom) panels show the samples binned in ipg (b/a). 
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Figure 12. The difference between the shape of the attenuation curve de- 
rived in this paper for galaxies with log 10 ips= —9.5 and b/a= 0.6, and 
the COO (stars) and MW (crosses) curves. The top (bottom) panel shows the 
comparison with the attenuation curve for low (high) /i* galaxies. The NUV 
point on the MW curve has been convolved with the GALEX NUV filter for 
this comparison, this convolution is not necessary for other wavelengths or 
for the COO curve (see Figure[8}. 



and ips according to simple linear functions: 



si = ci + c 2 (b/a) c + c 3 (ip c ) (22) 

52 = c 4 + c 5 (b/a) c + c 6 (ip c ) (23) 

53 = or + c 8 (b/a) c + cg(ip c ) (24) 
s 4 = 1.6 (25) 



where (b/a) c — (b/a) — 0.6 and ip c = ips + 9.5. The func- 
tions are valid in the primary parameter ranges covered by the 
binned samples studied in this paper: 0.3 < b/a < 0.9 and 
-10.2 < tp s < -9.3 for high fi* and -10.0 < ips < -9.1 fol- 
low fi* , and should not be extrapolated beyond these values. For 
galaxies outside of these ranges the functions could be calculated 
at the minimal or maximal extent of the range probed in our sam- 
ples. The free parameters for high and low fi* galaxies are provided 
in Tabled 

Figure [TO] shows the shape of the attenuation curves (Q\ ) for 
different values of ips and axis ratio. Figure [TJJ compares the ob- 
served attenuation in our binned galaxy samples to the appropriate 
empirical curves for the median ips and axis ratio of each sample. 
Overall the agreement is excellent, with no trend in residuals with 
either parameter. 

6.1 Comparison with other curves 

In Fig.[l2]we compare the dust curves derived here, to the MW and 
COO dust curves, for galaxies with ips~ —9.5 and b/a= 0.6. For 
a galaxy with attenuation in the stellar continua of Av = 1, the 
y-axis represents the difference in attenuation in magnitudes. The 
very steep attenuation in the blue found for the high /i* galaxies is 
evident here, as is the strong 2175A feature compared to the COO 
curve. For the low fi* galaxies, the attenuation is almost indistin- 
guishable to COO redwards of 5500A, although the curves differ 
substantially in the blue, and the lack of a 2 175 A feature in the 
COO curve is again evident. However, it is clear that for the low /i* 
galaxies, the feature is not as strong as in the MW. 



A complete dust attenuation curve is built from the shape 
(Q\), and the normalisation i.e. the amount of dust attenuation suf- 
fered by the continuum at 5500A (Tvr°l The most appropriate es- 
timator for tv will depend on the particular dataset. If the Balmer 
emission lines are used, Eqns. [13] and [J_5] provide the conversion 
between Ty t u ne and Tv, CO nt when the Balmer emission lines and 
continuum measurements are made within the same aperture. This 
is the best estimate possible with the SDSS data, but it assumes that 
the ratio measured within a 3 "diameter fibre is representative of the 
entire galaxy. Other possible estimations of ry c ould come from 
the m easurement of the dust emission in the FIR ( Ida Cunha et al. 
2008), or depth of the mid-IR silicate absorption feature Jwild et all 

bonh . 



7 DISCUSSION AND CONCLUSIONS 

Our main results are: 

NIR slope: The slope of the attenuation curve in the NIR is con- 
sistent with that measured for the MW extinction curve (~ 1.6) 
for all types of galaxies, irrespective of ips, inclination angle or the 
presence or absence of a bulge. 

Optical/UV slope: We find significant trends in the optical/UV 
slope of the attenuation curve with galaxy properties. The high /i* 
galaxies show a steeper blue/UV slope than the low galaxies 
and the average MW extinction curve. A flattening of the slope 
with increasing inclination is evident in both the high and low ^t* 
galaxies. There is no trend in optical slope with total optical depth. 

2175A feature: Most samples of galaxies exhibit a discontinu- 
ity in the UV attenuation curve that is suggestive of a 2 175 A dust 
feature at a level that is slightly weaker than seen in the MW. In the 
high n* galaxies the strength of the feature varies with axis ratio, 
from as strong as that observed in the MW extinction curve at low 
axis ratios, to undetected at high axis ratio. There is some evidence 
that the bump weakens with increasing ips- 

UV dust attenuation: For low (high) /1* galaxies with Av = 1, 
we find typically 0.3 ( 1 .0) magnitudes more attenuation in the NUV 
compared to the COO curve for starburst galaxies. 

Emission lines: We find that Balmer emission lines experience 
2-4 times more attenuation than the continuum at 5500A within 
the 3" SDSS fibre aperture. The ratio of emission line to contin- 
uum dust optical depths (Tvjme/Tv^ont) varies strongly with galaxy 
properties. It decreases with increasing ips, increases with increas- 
ing axis ratio (decreasing inclination) and is lower in high-/i* than 
in low-^i* galaxies. 

Radial gradients: Radial gradients are significant in all sam- 
ples, and increase with increasing ips - The strength of the gradients 
is not modified by the presence of a bulge. 



7.1 Comparison to the two component dust model 

The primary aim of this paper is to measure the shape of the dust 
attenuation curve of star-forming galaxies of different types. We 
now provide simple qualitative arguments to interpret the observed 



1 Note that this is different from Calzetti et al. who define the normali- 
sation by the amount of dust in the Balmer emission lines, and include the 
conversion between the amount of dust in the lines and continuum in the 
formula for the attenuation curve itself. 
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Figure 13. Pictorial representation of the distribution of dust in galaxies, which qualitatively can account for most results found in this paper. The bulge has 
been omitted for clarity, but is expected to contain predominantly old stars with little diffuse dust. 



trends in the context of the distribution of dust and stars in galax- 
ies. Fig.ll3lshows a pictorial representation of one possible config- 
uration for the global distribution of dust in star-forming galaxies, 
which can account for the majority of our results. This model for 
the distribution of dust relative to stars is based upon the popu- 
lar two component dust model, in which star-forming regions are 
surrounded by dense birth clouds, which evaporate on timescales 
similar to the main sequence lifetimes of B stars (~ 10 7 years). 
The remainder of the interstellar medium contains diffuse dust. 

This two component dust distribution h as been included 
in several models for dust in galaxies ( Silva et al 1 19981 ; 
ICharlot & Fallll2000l : |Popescu et alj|2000l ; lTuffs et alj|2004) . How- 
ever, many important details of the model are unconstrained by ob- 
servations. For example, how is the diffuse dust distributed relative 
to stars in the galaxy disks, how does dust content change with star 
formation rates, and what is the covering factor of the stellar birth 
clouds? Here, we attempt to address these questions in light of our 
new results combined with results in the literature. 



well with those estimated from mid-IR and far-IR radiation; nor by 
the weak trend between ips and b/a mentioned in Section[4] nor by 
the wea k trend between Ha equivalent width and inclination mea- 
sured bv lYip et al.l ( 120101) . Additional evidence against a significant 
skin effect comes from measurements of Ha rotation curves, which 
find less than uni ty optical depths to Ha at all but the steepest incli- 
nation angles jGoad & Robertslll98ll ; lGiovanelli & Havnesll2002l : 
lMisiriotisil2005l) . 

These observational results suggest that optically thick birth 
clouds have a low covering factor in ordinary star-forming galax- 
ies. Only in highly inclined galaxies do we see some evidence for a 
skin effect, as we observe a slight decrease in observed radial gra- 
dient strength (Figure[6). It is well known that optical data alone is 
insufficient to study some unusual galaxies, due to large dust opti- 
cal depths. Does this limitation of optical data apply to the general 
galaxy population? Our results suggest that only in highly inclined 
galaxies and galaxies with very high ips, where the diffuse dust op- 
tical depths are high due to long pathlengths or large diffuse dust 
contents, is optical data limited by dust. 



7.1.1 Covering factor of birth clouds 

As discussed in detail in SectionQ] one popular method for study- 
ing dust attenuation in galaxies, is by comparing the SEDs of in- 
clined and face-on galaxies. Through such a studv. rYTp et al.l J2010h 
found that the distribution of the ratio of Ha to H/3 luminosity does 
not vary with galaxy inclination. There are two plausible explana- 
tions for this result: (i) Emission lines are attenuated predominantly 
by the birth cloud from which they originate. The increase in dust 
opacity caused by an increase in path length through the diffuse 
ISM dust is insignificant in comparison to the difference in dust 
opacity between the birth clouds and ISM. (ii) Sightlines to star- 
forming regions that lie deeper within a galaxy are entirely opaque, 
causing us to measure the dust opacity only towards star-forming 
regions in an outer "skin". 

These different scenarios can alternatively be expressed in 
terms of the covering factor of stellar birth clouds in galaxies. If 
the covering factor is high, optical radiation from stars and gas pro- 
duced in the inner layers of the galaxy will have near-unit prob- 
ability of being absorbed by a dense cloud before escaping from 
the galaxy (skin model). Such an effect is not supported in general 
by the fact that SFR estimates from optical emission lines correlate 



7.1.2 Dust scale height and mixing with stars 

If stellar birth clouds have a small covering factor, most of the 
light dominating the UV-NIR emission from a star-forming galaxy 
has been attenuated by diffuse ISM dust alone. The geometrical 
arrangement of stars mixed with a dusty ISM can lead to an attenu- 
ation curve that is flatter (i.e. greyer) than the underlying extinction 
curve. This is because, at any wavelength, the emergent radiation 
is dominated by photons that are emitted by sources lying behind 
dust columns with less than unity optical depth. Thus, emergent 
red light originates from deeper within the galaxy than blue light, 
and attenuation appears more constant with wavelength. In Sec- 
tion |5.1| and Figure|4]we found that, in contrast, galaxies with high 
fi* have an attenuation curve that is steeper than expected at UV- 
optical wavelengths. This could occur if older/redder stars suffer 
less attenuation on average than younger/bluer stars, for example if 
both the younger stars and diffuse dust have a smaller scale height 
than older stars. In the case of high n* galaxies, there could be an 
additional effect from the presence of a bulge: galaxy bulges tend 
to have older stellar populations than galaxy disks, and sightlines 
not obscured by the disk may have lower dust contents. 
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The observed universality of the slope of the NIR attenuation 
curve (Fig. [5}, with a value equal to that of the MW extinction 
curve, implies a universal size distribution for large dust grains, 
and additionally sets further constraints on the relative geometry 
of old stars and dust in galaxies. For the old stars to suffer screen- 
like dust extinction, they must predominantly lie outside of the dust 
disk. This could occur through heating/migration effects of stars in 
galaxy disks, or through the thick disk being built through minor 
mergers of satellite galaxies which cease star formation as soon as 
they loose their gas suppli es. 

IXilouris etail i 19991) and iDalcanton etakl (2004) have also 
suggested that diffuse interstellar dust is distributed throughout the 
disk with a smaller scaleheight than that of the stars, following ob- 
servations of a small number of local galaxies. A va riable scale 
height for dust is also included in the models of e.g. iTuffs et ail 
( 2004). Matching the trends presented in this paper to detailed mod- 
els of the dust and star distribution in galaxies will reveal quantita- 
tive values such as the relative scale heights of stars and dust and 
the total opacity in galactic disks, as a function of specific star for- 
mation rate. 

7.1.3 Diffuse and birth cloud dust 

In Section [5"l4l and Figure [7] we observe strong trends in the ratio 
of line-to-continuum dust opacity as a function of axis ratio, ips 
and /i* . As discussed above, within the framework of the two com- 
ponent dust model the greater attenuation observed in the nebular 
emission lines compared to the stellar continuum is attributed to the 
additional attenuation caused by the stellar birth clouds surround- 
ing the hot OB stars that produce the emission lines. The covering 
factor of these birth clouds is unknown, although we argued above 
that it is likely to be small. If this is the case then we can, to first 
approximation, assume that the attenuation of the emission lines 
arises primarily from dust in the birth clouds, and the attenuation 
of the stellar continuum arises from dust in the diffuse ISM. 

In more inclined galaxies, the pathlength through the dif- 
fuse dust is greater: the observed opacity towards all sources lying 
within the dust disk increases with increasing inclination. This can 
have a significant effect on the optical depth measured in the con- 
tinuum. By virtue of the significantly larger opacities of the birth 
clouds compared to the diffuse ISM dust, inclination will have less 
impact on the total optical depth measured in the Balmer emission 
lines so long as the opacity of the ISM is low enough. Thus, the 
observed correlation between the ratio of line-to-continuum dust 
opacity and axis ratio, and the inva riance of ratio of Ha to H/3 lu- 
minosity with inclination found by I Yip et al.l ( |2010|) . is consistent 
with a two-component dust model, in which the birth clouds have 
significantly greater dust optical depths than the diffuse ISM. The 
increase in strength of the 2 175 A absorption feature with inclina- 
tion in the bulge-dominated galaxies may also be related to an in- 
crease in diffuse dust opacity. However, the lack of a similar trend 
in the low /x* galaxies is puzzling. The answer to this question may 
provide clues to the nature of the grain responsible for the 2 175 A 
absorption feature. 

At least two factors could cause the trend found between the 
ratio of line-to-continuum dust opacity and ips. Firstly, the greater 
contribution of young stars to the V-band continuum light in high 
ips galaxies may cause the measured Tv.cont to approach Tv.iine at 
high ips if these younger stars suffer greater dust attenuation rela- 
tive to older stars. Secondly, a decrease in diffuse dust content from 
star-forming to quiescent galaxies would reproduce the observed 
trend. Given the significant amplitude of the observed trend, and 



the small quantities of diffuse dust found in elliptical galaxies, we 
believe this to be the most plausible explanation. Such an evolution 
would have implications for theories of dust grain formation dur- 
ing star formation, and dust grain destruction in the diffuse ISM. 
The larger ratio of line-to-continuum dust opacity in low fi* galax- 
ies may relate to a decreased diffuse dust content in these lower 
metallicity galaxies. The trend is in qualitative agreement with ob- 
servations of diffuse FIR emission, which show that later type spi- 
rals having a lower fract ion of diffuse dust than early type spirals 
dSauvage & Thuanll 19921) . 

Finally, the strong radial gradients found in the attenuation of 
the stellar continua (Fig.[6j imply that the central regions of galax- 
ies are dustier than the outskirts, a result presumably related to the 
metallicity gradients observed in galaxy disks. The link between 
ips and dust content is intriguing, but observationally we are lim- 
ited by the 3" SDSS fibre; spatially resolved spectroscopic studies 
of spiral galaxies will allow us to explore this area further. 

7.1.4 Summary 

Taken together, our observations are broadly consistent with the 
popular two-component model for the distribution of dust in galax- 
ies, although several points remain unexplained. In Fig.[T3]we have 
suggested some improvements to the basic picture, based upon re- 
sults in this paper. In detail, we suggest a low covering factor for 
the dense birth clouds surrounding the youngest stars, we infer a 
scale height for the dust that is smaller than that of the old stars, 
and we suggest a radial gradient in the amount of diffuse ISM dust 
with more dust found towards the center of the galaxy, particularly 
in galaxies with high ips. Clearly, a more quantitative comparison 
between our results and detailed models will lead to a deeper under- 
standing of the production and destruction of dust in the Universe. 

7.2 Applying the dust attenuation curve 

The general properties of the galaxies studied in this paper are sum- 
marised in Figure [TJ and the dust curves presented here should be 
applied with caution to galaxies outside of these ranges. In order to 
apply the dust attenuation curves provided in this paper to the SED 
of a star-forming galaxy: 

• Does the galaxy have a significant bulge? If so, use the coef- 
ficients for the high fi* galaxies. Otherwise use the coefficients for 
the low fj* galaxies. 

• Determine the inclination and ips of the galaxy, if either are 
unknown set b/a= 0.6 and logips — —9.5. 

• Determine the Tv, CO nt of the galaxy. If using the Balmer emis- 
sion lines: 

- Account for strong dependence of TvjineM/.cont on inclina- 
tion and ips. 

- Ensure that the emission lines are measured using the same 
aperture as the continuum data, or correct for radial gradients 
using Figure [6] 

• Determine the shape of the appropriate attenuation curve from 
Eqn.[T8] or using the IDL code provided in Appendix [C] 

The attenuation curves presented in this work, like that of 
COO, account for the effects of mixed dust-star geometry and global 
geometry in galaxies. Unlike the MW, or Magellanic extinction 
curves, our attenuation curve does not assume a screen-like dust 
geometry. Additionally, including the dependence on ips accounts 
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for the different forms of attenuation that different types of stars 
suffer, and the changing balance of these types of star to the inte- 
grated light of the galaxy depending on the SFR of the galaxy. 
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APPENDIX B: SQL QUERIES 

The queries below are run on the SDSS Catalog Archive 
Server (http://casjobs.sdss.org/casjobs/ 1, WFCAM Science Archive 
(http://surveys.roe.ac.Uk/wsa/i and MAST Catalog Archive Server 
(http://mastweb.stsci.edu/gcasjobs/! respectively. The resulting 
data tables are joined manually on SDSS specobjid, using a pur- 
pose designed IDL code. Here we provide an extract of each query, 
the full SQL code is available online. 



APPENDIX A: TESTS OF THE PAIR-MATCHING 
METHODOLOGY 

As described in Section l2~2l there are some important features of the 
new pair-matching method to be aware of if the same experiment 
is to be repeated on different datasets. In Figures [ATI and lA2l we 
present some simple tests to show that our method is valid for the 
dataset in this paper. 
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Figure Al. The optical slope (left) and Axvjine^fv.cont ra ti° measured from samples of low (top panel) and high (bottom panel) fi* galaxies, where the 
pairs of galaxies are required to have ATv,ii n e= [0.3, 0.5, 0.7, 0.9] ± 0.2 (points from left to right). This shows that there is no variation in the shape of the 
dust curves or in the line-to-continuum optical depth ratio, with continuum optical depth. 




Figure A2. On the left we compare attenuation curves derived from galaxy pairs matched on ips estimated from Ha within the 3"SDSS fibre divided by 
stellar mass estimated from the optical photometry (left), and a pure photometric estimate of ips (right). We find that the observed attenuation curves are 
substantially weaker when using the pair matching based upon photometric estimates of ipg. On the right we investigate the reason for this, comparing the 
two SFR estimates as a function of ratio of Ha to H/3 luminosity (2D histogram in log number of galaxies). We find that the photometric SFR= M*ipgis 
systematically lower in galaxies with a high ratio of Ho to H/3 luminosity. This is not unexpected as it is difficult to distinguish dust reddening from an older 
stellar population using optical photometry alone. The resulting correlation between stellar population and dust content causes the attenuation curve to appear 
considerably flatter. 
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/* To obtain parameters for spectroscopic galaxies in the SDSS */ 

/* Basic parameters */ 

SELECT s.bestobjid,s. specob j id, m jd, plate , fiberid, s . z , sn_0 , sn_l , 
petror50_r , petror90_r, petror50_z , rowindex, p . extinction_g, p . petrorad_r , 

/* Petrosian magnitudes */ 

p . petromag_u, p . petromag_g, p . petromag_r, p . pet romag_i , p . pet romag_z , 

/* The number of annuli for which there is a measurable signal */ 
p.nprof_u, p.nprof_g, p.nprof_r, p.nprof_i, p.nprof_z, 

/* The u-band AperFlux */ 

upO.profMean as uAperFluxO, upl.profMean as uAperFluxl, 

up2.profMean as uAperFlux2, up3.profMean as uAperFlux3, 

up4.profMean as uAperFlux4, up5.profMean as uAperFlux5, 

up6.profMean as uAperFlux6, up7.profMean as uAperFlux7, 

/* The g-band AperFlux */ 

gpO.profMean as gAperFluxO, gpl.profMean as gAperFluxl, [....] 

/* mydb . sdssmpa_rowind contains the SDSS-MPA rowindex and specob jid for all galaxies */ 
FROM specob j as s, photoobjall as p, mydb . sdssmpa_rowind as m, 

PhotoProfile as upO, PhotoProfile as gpO, PhotoProfile as rpO, PhotoProfile as ipO, 
PhotoProfile as zpO, PhotoProfile as upl, [....] 



/* Join Specob j, PhotoobjAll, SDSSMPA_rowind and select only spectroscopic galaxies */ 
WHERE s . specob jid=p . specob jid AND s . specob jid=m. specob jid AND s . specclass=2 AND 



/* Join u-band PhotoProf iles */ 

s.bestobjid = upO.objid AND upO 

s.bestobjid = upl.objid AND upl 

s.bestobjid = up2.objid AND up2 

s.bestobjid = up3.objid AND up3 

s.bestobjid = up4.objid AND up4 

s.bestobjid = up5.objid AND up5 

s.bestobjid = up6.objid AND up6 

s.bestobjid = up7.objid AND up7 



.band = AND upO.bin = AND 

.band = AND upl. bin = 1 AND 

.band = AND up2.bin = 2 AND 

.band = AND up3.bin = 3 AND 

.band = AND up4.bin = 4 AND 

.band = AND up5.bin = 5 AND 

.band = AND up6.bin = 6 AND 

.band = AND up7.bin = 7 AND 



/* Join g-band PhotoProf iles */ 
s.bestobjid = gpO.objid AND gpO.band 



1 AND gpO.bin = AND 
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/*To obtain parameters for UKIDSS sources which are spectroscopic galaxies in the SDSS. */ 

SELECT ss . specob jid, ls.sourceid, x . distancemins, 

/* Petrosian magnitudes */ 

Is . yPetroMag, Is . j_lPetroMag, Is . hPetroMag, Is . kPetroMag, 



/* dust attenuation */ 
Is.aY, ls.aH, Is.aJ, ls.aK, 



/* Petrosian radii */ 

yd.PetroRad as yPetroRad, jd.PetroRad as j_lPetroRad, hd.PetroRad as hPetroRad, 
kd.PetroRad as kPetroRad, 

/* UKIDSS Aper mags */ 

yd.Aperfluxl as yAperFluxl, jd.AperFluxl as j_lAperFluxl , hd.AperFluxl as hAperFluxl, 
kd.AperFluxl as kAperFluxl, 
yd.Aperf lux2 as yAperFlux2, [....] 



/* J band exposure time */ 
mf . expTime 



FROM lasYJHKSource AS Is, bestdr7 .. specob j as ss, lasSourceXDR7Photo0b j as x, 
lasY JHKMergeLog AS ml, lasDetection as jd, lasdetection as hd, lasdetection as yd, 
lasdetection as kd, multiframe as mf 

WHERE x .masterob jid=ls . sourceid AND ss . bestob j id=x . slaveob j id AND 
Is . f ramesetID = ml . frameset ID AND ml.j_lmfID = mf .mult if ramelD AND 

/* Join merge log to detection for the J frame: */ 

ml.j_lmfID = jd.multif ramelD AND ml . j_leNum = jd.extNum AND 

/* Join source to detection to look up the detection required: */ 

ls.j_lSeqNum = jd.seqNum AND 

/* Join merge log to detection for the Y frame: */ 

ml.ymfID = yd.multif ramelD AND ml . yeNum = yd.extNum AND 

/* Join source to detection to look up the detection required: */ 

ls.ySeqNum = yd.seqNum AND 



/* Join merge log to detection 
ml.hmfID = hd.multif ramelD AND 
/* Join source to detection to 
ls.hSeqNum = hd.seqNum AND 

/* Join merge log to detection 
ml.kmfID = kd.multif ramelD AND 
/* Join source to detection to 
ls.kSeqNum = kd.seqNum AND 



for the H frame: */ 

ml . heNum = hd.extNum AND 

look up the detection required: */ 

for the K frame: */ 

ml.keNum = kd.extNum AND 

look up the detection required: */ 



/* remove blended objects */ 

Is . yppErrBits&OxOOOOOOlO = AND Is . j_lppErrBits&0x00000010 = AND 
Is . hppErrBits&OxOOOOOOlO = AND Is . kppErrBits&OxOOOOOOlO = AND 

/* now get spectroscopic galaxies and select the nearest match */ 
ss . specclass=2 AND distanceMins<0 . 033333 AND 

distanceMins IN ( SELECT MIN (distanceMins ) FROM lasSourceXDR7PhotoOb j 
WHERE masterObjID=x.masterObjID AND sdssPrimary=l ) 
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/*To obtain parameters for GALEX sources which are spectroscopic galaxies in the SDSS. */ 

/* Basic parameters */ 

SELECT ss . bestob j id, ss . specob jid, gg.objid, gg. distance, gg . reversemultiplematchcount , 
/* GALEX Fluxes */ 

ga.nuv_flux, ga . nuv_f luxerr , ga . f uv_f lux, ga . f uv_f luxerr , ga . e_bv 

/* myDB . dr7SFgal contains an uploaded list of all the specob jids we want to match */ 
FROM myDB . dr7SFgal as ss, galexgr4plus5 . . xsdssdr7 as gg, galexgr4plus5 . .photoob j as ga 

WHERE ss . bestob jid=gg . SDSSobj id and gg . ob j id=ga . ob j id and gg . reversedistancerank=l 
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FUNCTION VW_DUSTCURVE, wave, flux, tau_v, flag_mustar, ssf r=ssf r, ba=ba, 

silent =silent, magnitudes=magnitudes 

+ 

NAME: VW_DUSTCURVE 

AUTHORS : 

Vivienne Wild <vw@roe . ac . uk> 

PURPOSE : 

Calculate the unattenuted galaxy spectral energy distribution 
followingt the prescription of Wild et al . 2011 

CALLING SEQUENCE: 

unattenuated_f lux = VW_DUSTCURVE (wave, flux, tau_v, f lag_mustar , [ssfr=, ba=] ) 
INPUTS : 

wave = wavelength array in microns, not to extend beyond [0.14,2.05] 

flux = flux array (or magnitudes, see /magnitudes keyword) 

tau_v = optical depth in the stellar continuum at 5500Angstrom 

flag_mustar = -1 for high stellar surface mass density, or bulge-dominated 
= 1 for low stellar surface mass density, or disk-dominated 

OPTIONAL INPUTS: 

ssfr = logarithm of the specific star formation rate of 

the galaxy in years ( log (psi_s/yr " -1 ) . If none is 
supplied, -9.5 is assumed 
ba = minor/major axis ratio of galaxy. If none is 

supplied then 0.6 is assumed. 

KEYWORD PARAMETERS : 

silent = set to prevent some messages 

magnitudes = set to correct an array of magnitudes, rather than fluxes 
OUTPUTS : 

unattenuated_f lux = unattenuated flux array 
REFERENCES: [1] Wild, Chariot, Brinchmann et al . , 2011, MNRAS, submitted 
NOTES : 

MODIFICATION HISTORY: 

2011 First implementation in IDL V. Wild 

********************************************* 
Copyright (c) 2011, Vivienne Wild 

Permission to use, copy, modify, and/or distribute this software for any 
purpose with or without fee is hereby granted, provided that the above 
copyright notice and this permission notice appear in all copies. 

THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES 
WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF 
MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR 
ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES 
WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN 
ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF 
OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. 
****************************************************************************************; 
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; Return to caller on error. 
On_Error, 2 

n = 20d ; smoothing parameter 

lc_l_eff = 0.2175d ; rest-frame break points 

lc_2_eff = 0.3d 

lc_3_eff = 0.8d 

V = 0.55 ; normalised here 

;;— verify inputs 

if N_params() LT 4 then begin 

print, 'syntax: VW_DUSTCURVE (wave, flux, tau_v, flag_mustar, [ssfr=, ba=] ) ' 

return, -1 
endif 

;;— check for b/a and ssfr 
if n_elements (ssf r) eq then begin 
ssfr=-9.5 

if not (keyword_set ( silent ) ) then print, 'Assuming log (psi_s ) =-9 . 5 ' 
endif 

if n_elements (ba) eq then begin 
ba = 0.6 

if not (keyword_set ( silent ) ) then print, 'Assuming b/a=0.6' 
endif 

; ; — define parameters for high / low mustar cases 

if flag_mustar eq -1 then P = [1.3, 0.6, 0.1, 1.1, 0.3, -0.2, 0.2, 0.9,0.1] $ 
else if flag_mustar eq 1 then P = [1.1, 0.4, -0.1, 0.7, 0., 0.4, 0.15, 0.,0.2] $ 
else begin 

print, 'Please specify: flag_mustar = -1 for high stellar surface mass density' 
print, ' = 1 for low stellar surface mass density' 

print, 'syntax: VW_DUSTCURVE (wave, flux, tau_v, flag_mustar, [ssfr=, ba=] ) ' 
return, -1 
endelse 

; ; — check parameters are within specified limits 
ind = where (wave It 0.14 or wave gt 2.05) 
if ind[0] ne -1 then begin 

print, 'Dust curve only valid in wavelength range . 14-2 . 05microns' 

return, -1 
endif 

if ba It 0.3 or ba It 0.9 then begin 

print, 'Dust curve only valid in axis ratio range 0.3-0.9, see Section 6 of [1]' 

return, -1 
endif 

if flag_mustar eq -1 then begin ;high mustar 
if ssfr It -10.2 or ssfr It -9.3 then begin 

print, 'For bulge-dominated galaxies dust curve only valid in ssfr range -10.2 -> -9.3' 
return, -1 
endif 
endif 

if flag_mustar eq 1 then begin ; low mustar 

if ssfr It -10.0 or ssfr It -9.1 then begin 

print, 'For disk-dominated galaxies dust curve only valid in ssfr range -10.0 -> -9.1' 
return, -1 
endif 
endif 
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; ; — define relations with SSFR and b/a 

ba_c = ba - . 6 
ssfr_c = ssfr + 9.5 

; ; sopt 

sopt = P [0] +ba_c*P [1] +ssfr_c*P [2] 

; ; snir 
snir = 1.6 

; ; snuv 

snuv = P [3] +ba_c*P [4] +ssfr_c*P [5] 
; ; sfuv 

sfuv = P [ 6] +ba_c*P [7] +ssf r_c*P [8] 



lc_l = lc_l_eff 

lc_2 = ( (lc_l " snuv) / (lc_2_ef f " (snuv-sopt ) ) ) " (1/sopt ) 
lc_3 = ( (lc_2 " sopt ) / (lc_3_ef f " (sopt-snir ) ) ) " (1/snir) 

norm = ( ( (V/lc_l) " (sfuv*n) + (V/lc_l) " (snuv*n) + (V/lc_2) " (sopt*n) +$ 
(V/lc_3) " (snir*n) ) " (-1/n) ) ; Qlambda=l at V 

Q_lambda = ( ( (wave/lc_l ) " ( sf uv*n) + (wave/ lc_l ) " ( snuv*n) + (wave/lc_2 ) " (sopt *n) +$ 
(wave/lc_3 ) " ( snir *n) ) " (-1 /n) ) / norm ;eq. 18 

tau_lambda = Q_lambda*tau_V 

if not (keyword_set (magnitudes) ) then return, f lux/exp (-tau_lambda) $ ;e.g. eq 5 
else return, flux - 1 . 086*tau_lambda ; eqn . 3 

END 
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